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ABSTRACT 

We present optical photometry and spectroscopy of five type la supernovae discovered by the Nearby Su- 
pernova Factory selected to be spectroscopic analogues of the candidate super-Chandrasekhar-mass events 
SN 2003fg and SN 2007if. Their spectra are characterized by hot, highly ionized photospheres near maximum 
light, for which SN 1991T supplies the best phase coverage among available close spectral templates. Like 
SN 2007 if, these supernovae are overluminous (-19.5 < My < -20) and the velocity of the Si II A6355 ab- 
sorption minimum is consistent with being constant in time from phases as early as a week before, and up to 
two weeks after, Z?-band maximum light. We interpret the velocity plateaus as evidence for a reverse-shock 
shell in the ejecta formed by interaction at early times with a compact envelope of surrounding material, as 
might be expected for SNe resulting from the mergers of two white dwarfs. We use the bolometric light curves 
and line velocity evolution of these SNe to estimate important parameters of the progenitor systems, including 
""^Ni mass, total progenitor mass, and masses of shells and surrounding carbon/oxygen envelopes. We find 
that the reconstructed total progenitor mass distribution of the events (including SN 2007if) is bounded from 
below by the Chandrasekhar mass, with SN 2007if being the most massive. We discuss the relationship of 
these events to the emerging class of super-Chandrasekhar-mass SNe la, estimate the relative rates, compare 
the mass distribution to that expected for double-degenerate SN la progenitors from population synthesis, and 
consider implications for future cosmological Hubble diagrams. 

Subject headings: white dwarfs; supernovae: general; supernovae: individual (SN 2003fg, SN 2007if, 
SN 2009dc, SNF 20080723-012) 
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Type la supernovae (SNe la) have become indispensable as 
luminosity distance indic ators for exploring the accelerated 
expan sion of the universe (iRiess et a/.lll998l ; lPerImutter et al.\ 
II999I) . Their utility is due mainly to their very high lumi- 
nosities, in combination with a set of relations between in- 
trinsic luminosity, color and light curve width ("Riess et o/J 
II998; Tripp 1998; Phillips et al. 1999; Goldhaber et al. 200 J 
which reduces their dispersion around the Hubble diagram to 
^ 0.15 mag. Much recent attention has been given to im- 
proving the precision of distance measurements by search- 
ing for further standardization relations, with some meth- 
ods using near-maximum-light spectra to deliver core Hubble 
residual dispersions as low as 0.12 mag (Bailev et al. 2009i 
IWang et al. 2009; Folatelli et al. 2009; Foley & Kasen 201 ij. 

Despite ongoing research, however, many uncertainties re- 
main regarding the physical nature of SN la progenitor sys- 
tems, althoug h observational subclasses can be formed (e.g., 
iBranch et a?TlI993t iBenetti et al\ 120051) . Detailed observa- 
tions of the light curve of the nearby type la SN 201 Ife shortly 
after explosion have shown that it must have had a compac t 
progenitor (PTFl Iklv lNugent et al\20\ UlBloom et a/.l20I2l) . 
in line with expectations that SNe la result from thermonu- 
clear explosions of white dwarfs in binary systems. The mass 
and evolutionary state of the progenitor's companion star, the 
events that trigger the explosion, and the circumstellar and 
host galaxy environment of most normal SNe la remain un- 
known. Next-generation SN la cosmology experiments make 
stringent demands, and any cosmological or astrophysical 
phenomenon which could bias the measured luminosities of 
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SNe la at the level of a few p ercent has become important 
to mvestigate ( lKimefa/.ll200l . If two or more progenitor 
channels exist corresponding to different peak luminosities or 
luminosity standardization relations, evolution with redshift 
of the relative rates of observed SNe la from these channels 
could mimic the effects of a time-varying dark energy equa- 
tion of state (Under 2006). 

The two main competing SN la progenitor scenario s are 
the single-degenerate scenario (IWhelan & IbenI 119731) . in 
which a carbon/oxygen white dwarf slowly accretes mass 
from a non-degenerate companion until exploding near the 
Chandrasekhar mass, and the double-degenerate scenario 
dlben & Tutukovl 119841) . in which two white dwarfs col- 
lide or merge. There are also sub-Chandrasekhar models 
jWooslev & Weaved [19941: Ivan Kerkwiik. Chang. & JusthamI 
120101) . in which the explosion is triggered by the detonation 
of a helium layer on the surface o f a sub-Chandrasekhar- 
mass white dwarf jSim et al\ 120101) . Although historically 
the single-degnerate scenario has been favored, the double- 
degenerate scenario has rec ently gained ground both th eoret- 
ically and observationally. iGilfanov & BogdanI (120101) used 
X-ray observations of nearby elliptical galaxies to set a strin- 
gent limit on the number of accreti ng white dwarf systems 
in those galaxies (see also lDi Stefanoll20103) . and hence on 
the single-degenerate contribution to the SN la rate, although 
their interpretation of the measurements has been questioned 
jDi Stefanoll2010bl: iHachisu et a/.ll2010l) . Based on searches 
for an ex-companion star near the site of explosion, iLi et al\ 
(1201 Ih ruled out a red giant companion for single-degenerate 
models of SN 201 Ife, and Schaefer & Pagnotta (2012) argue 
strongly that the supernova remnant SNR 0509-67.5, in the 
Large Magellanic Cloud, must have had a double-degenerate 
progenitor; however, single-degenerate scenarios have been 
put forth in which long delays between formation of the pri- 
mary white dwarf and the SN la explosion could allow the 
companion to evolve a nd become fainter, evadin g attempts 
to detect them directly dPi Stefano & Kilicll20T2h . While a 
merging white dwarf system could also undergo accretion- 
induced collapse to a neutron star dNomoto & Kondolll991l : 
ISaio & NomotolfT998h rather than exploding as a SN la, the- 
oretical investigations of SNe la from mergers have also pro- 
gressed. Some merger simulations produce too much unburnt 
materi al to reproduce spectra of normal SNe la dPfannes et al\ 
I2010al) : other simulations suggest that if the merger process 
is violent enough to ignite the white dwarf promptly, me rg- 
ers may produce subluminous SNe la (P akmor et a/.ll20ll1) or 
even normal SNe la (Pakmor et al. 2012). 

Interest has also been aroused by the discovery of "super- 
Chandra" S Ne la which far exc eed the norm i n luminosity: 
SN 2 003fg dHowellef a/.ll2006h. SN 20 06gz dHicken et al.\ 



dYan 



|2007|) . SN 2007if jScalzo ef aLl l20Tol: lYuan et al. 



20m 

and SN 2009dc dYamanaka ef a/. 2009; T anaka,eLa/J |2010|; 
iTaubenberger et al.WlOl iTiSilverman et a/.il201 lb . If the light 
curves of these SNe are powered by the radioactive decay 
of ^^Ni, as expected for normal SNe la, the ^^Ni mass nec- 
essary to produce the observed luminosity implies a sys- 
tem mass significantly in excess of the Chandr asekhar mass. 
The high luminosity of SN 2006gz reported in i Hicken et al.\ 
d2007i) hinges upon an uncertain reddening correction, and 
late-phase photometry and spectroscopy suggest a smaller 
^^Ni mass tha n that inferred from the dereddened peak lu- 
minosity (Ma edaef a/.ll2009h : SN 2003fg, SN 2007if, and 
SN 2009dc are much more luminous than normal SNe la even 



before dereddening. None of these SNe lie on the existing lu- 
minosity standardization relations. If clear photometric and 
spectroscopic signatures of the progenitor system or explo- 
sion mechanism can be discovered for these rare SNe la, they 
may help bring to light similar, but weaker, signatures in less- 
extreme SNe la. 

In addition to its high luminosity (My = -20.4), SN 2007if 
showed a red {B-V = 0.18) color at maximum light, C II 
absorption in spectra taken near maximum light, and a low 
(~ 8500 km s~') and very slowly-e volving Si II A6355 ab- 
sorption velocity. IScalzo et al.\ d2OI0l) used this information to 
model SN 2007if as a "tamped detonation" resulting from the 
explosion of a super-Chandrasekhar-mass white dwarf inside 
a dense, but compact, carbon/oxygen envel ope, as expected 
in some double-degenerate merger scenarios dkhokhlov et al.\ 
[1993; Hoflich & Khohklov 1996). Using the bolometric light 
curve between 60 and 120 days past explosion to determine 
the optical depth for gamma-ray trappi n g in the ejecta (see 
e.g. lJeffervlll"999l:lStritzinger et a/.ll2006l) . (Scalzo et aIl(l207ol) 
estimated the total SN 2007if system mass to be 2.4 ± 0.2 Mq, 
with ^ 15% of this mass bound up in the carbon/oxygen en- 
velope formed in the merger process. Such a high mass is 
near the theoretical upper mass limit for two carbon/oxygen 
white dwarfs in a binary system. However, numerical calcu- 
lations have since confirmed that the very large mass of ^''Ni 
needed to explain the luminosity of SN 2007if can plausi- 
bly be produced in a collision of two high-mass white dwarfs 
([Raskin et al. 2010), or in the prompt detonation o f a single 
rapidly -rotating white dwarf (iPfannes et al.ll2010b[) . Recent 
work has suggested that a white dwarf with a non-degenerate 
companion could be spun up by accretion to high mass, and 
remain rotationally suppo rted for some t ime only to explode 
later (Justham 201 1; Hac hisu et a/.ll201I[) . 

In contrast, the comparably-bright SN 2009dc had a nor- 
mal color (B - V = 0) near maximum light and showed 
rapid Si II velocity e volution at ~ 100 km s"' day"' 
dSilverman ef g/.l 120111). d i fficult to explain by a tamped 
detonation. iTanaka et al] (120101) group SN 2009dc with 
the 2003fg-like SNe la based on its high luminosity and 
broad light curve, but they compare it spectroscopically with 
SN 2006gz, given its stro ng Si II A6355 and C I I A6580 ab- 
sorption at early phases. [Silverman et al.\ (1201 ll) also noted 
some spectroscopic differences in the post-maximum spectra 
of SN 2007if and SN 2009dc. Moreover, the SN was fainter at 
one ye ar after explosion than expected for the estimated^ ^Ni 
mass d Taubenberger et al.\ 1201 U [Silverman et al.' '201 11), as 
iMaedaef g/. (2009) noted for SN 2006gz. .Taubenbergeref al.\ 
d2011l) calculated a total system mass of 2.8 Mq for 
SN 2009dc by estimating the diffusion tim e from the width 
of the bolometric light curve dArnettlll982b : they explored a 
number of different thermonuclear and core-collapse explo- 
sion scenarios, and found none of them to be completely sat- 
isfactory in explaining all the observations. 

After the discovery of SN2007if, we remained vigilant 
for SNe with similar characteristics; when possible, addi- 
tional follow-up was obtained when such SNe were recog- 
nized. Given the high ionization state and predominance of 
Fe II and Fe III absorption in SN 2007if's spectra up until 
maximum light, characteristics shared with SN 199 IT, we 
used a 1991T-like spectroscopic classification to select for 
super-Chandrasekhar-mass SN candidates, triggering follow- 
up even for more distant, fainter examples. The Nearby Su- 
pernova Factory (SNfactory) obtained, as part of its spectro- 
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scopic follow-up program on a large sample of nearby SNe la, 
observations of five such SNe la besides SN 2007if, the pre- 
sentation of which is the subject of this paper. Later exam- 
ination of the full data set of SNfactory spectroscopic time 
series showed that the SNe in our sample each also show a 
plateau in the time evolution of the velocity of the Si II A6355 
absorption minimum, lasting from the earliest phase the ve- 
locity was measurable until 10-15 days after maximum light, 
as in SN 2007if. The absorption minimum velocities of other 
intermediate-mass elements in these SNe also show plateau 
behavior These events have a different appearance from 
SN 2006gz and SN 2009dc, which do not show velocity 
plateaus and show somewhat different behavior in their early 
spectra and late-time light curves. SN 2003fg was spectro- 
scopically observed only once, making it impossible to deter- 
mine how it evolved spectroscopically. 

Our supernova discoveries, our sample selection, and the 
provenance of our data are described in ^ the light curves 
and spectra are presented in section ^ In ^ we model 
our SNe as tampe d detonations, using the formalism of 
IScalzo et al\ (120101) to estimate an envelope mass and total 
system mass for each SN. We discuss the broader implica- 
tions of our results in |5j including implications for progen- 
itor systems, explosion mechanisms, and cosmology, and we 
summarize and conclude in ^ 

2. OBSERVATIONS 

This section details the discovery, selection and follow-up 
data for our sample of candidate super-Chandrasekhar-mass 
SNe la. 

2.1. Discovery 

The supernovae are among the 400 SNe la discovered in the 
SNfactory SN la search, car ried out between 20 05 and 2008 
with the QUEST-II camera ( iBaltav ef a/.ll2007h mounted on 
the Samuel Oschin 1.2-m Schmidt telescope at Palomar Ob- 
servatory ("Palomar/QUEST"). QUEST-II observations were 
taken in a broad RG-610 filter with appreciable transmission 
from 6100-10000 A, covering the Johnson R and / band- 
passes. Table[T]lists the details of the SN discoveries, includ- 
ing SN 2007if. 

Upon discovery candidate SNe were spectroscopically 
screened using the Sun erNo va Integral F ield Spectrograph 
rSNIFS: lAidering et amml iLantz ef a/.ll20()l on the Uni- 
versity of Hawaii (UH) 2.2 m on Mauna Kea. Our normal cri- 
teria for continuing spectrophotometric follow-up of SNe la 
with SNIFS were that the spectroscopic phase be at or be- 
fore maximum light, as estimat ed using a template- matching 
code similar e.g. to S UPERFIT [Howell ef a/.l(l2005h . and that 
the redshift be in the range 0.03 < z < 0.08. Aware of the 
potential for discovering nearby counte rparts to the candidat e 
super-Chandrasekhar-mass SN 2003fg iHowell et aT\ (Hool), 
we allowed exceptions to our nominal redshift limit for con- 
tinued follow-up if the spectrum of a newly-screened candi- 
date appeared unusual or especially early. 

2.2. Selection criteria 

Our spectroscopic selection was informed by the existing 
spectra of SN 2003fg and SN 2007if. In particular, the pre- 
maximum and near-maximum spectra of SN 2007if showed 
weak Si II A6355 and Ca II Hh-K absorption and strong absorp- 
tion from Fe II and Fe III, with noted similarity to SN 199 IT 
(IScalzo et a/.ll2010l) . We therefore prioritized spectroscopic 



follow-up for SNe la visually similar to, or more extreme (i.e. 
having weaker IMF and stronger Fe-peak absorption) than, 
SN 199 IT itself. Here we have excluded the less extreme 
1999aa-like s, which can be separat ed based on the strength of 
Ca II Hh-K (ISilverman et a/.ll2012l) . As a cross-check on our 
initial selection conducted usi ng the initial classificatio n spec- 
tra, we have run SNID v5.0 (iBlondin & Tonryll2007h on the 
entire SNfactory sa mple, using version 1.0 of the templates 
supplemented by the lScalzo et g/.l (120101) SNfactory spectra of 
SN 2007if. We searched for pre-maximum spectra for which 
the best subtype was "Ia-91T", or "la-pec" with the top match 
being a 1991T-like SN la or SN 2007if itself; this yielded the 
same set of SNe as the visual selection. 

Our sample could thus be described as 199 IT- like based on 
their spectroscopic properties. However, this categorization is 
often used to imply lightcurve characteristics — luminosity 
excess or slow decline rate — that were not part of our selec- 
tion criteria. Moreover, since our interest is ultimately in the 
masses of these systems, we will refer to these as "candidate 
super-Chandra SNe la" throughout this paper 

The sample of six objects (including SN 2007if) presented 
here constitutes all such spectroscopically-selected candidate 
super-Chandra SNe la in the SNfactory sample. This was es- 
tablished as part of the classification cross-check described 
above. An additional 141 SNe la discovered by the SNfactory 
were also followed spectrophotometrically and constitute a 
homogenous comparison sample. Like SN 2003fg, but unlike 
SN 2006gz and SN 2009dc, our sample of super-Chandra can- 
didates and our reference sample were discovered in a wide- 
area search and therefore sample the full range of host galaxy 
environments. This may prove important in understanding 
the formation of super -Chandra SNe la, e.g., if metallicity 
plays an important role ([Taubenberger et a/.l201 UlKhan et al\ 
1201 It iHachisu et al\ |20111k We reiterate that characteris- 
tics such as luminosity excess, velocity evolution, lightcurve 
shape, etc. were not used in our selection — it is purely based 
on optical spectrophotometry. 

2.3. Follow-up Observations 

Most of the BVRI photometry in this work was synthe- 
sized from SNIFS flux-ca l ibrated rest-frame spectra, using 
the bandpasses of iBesselll (119901) . an d corrected for Galac- 
tic dust extinction using E{B-V) fromlSchlegel erZl (ll998h 
and the extinction law of lCardelli et a/.l (11988!) with Rv = 3.1. 
Follow-up BVRI photometry for SNF 20070803-005, using 
the ANDICAM imager on the CTIO 1.3-m, was obtained 
through the Small and Moderate Aperture Research Telescope 
System (SMARTS) Consortium. Redshifts were obtained 
from host galaxy spectra, which were either extracted from 
the SNIFS datacu bes, or taken s eparately using the Kast Dou- 
ble Spectrograph (lMilleiill993l) on the Shane 3 m telescope 
at Lic k Observatory, t he Low Resolution Imaging Spectro- 
graph (lOkeef a/.ll995h at Keck-I on Mauna Kea or the Good- 
man High-Throughput Spectrograph at SOAR on Cerro Pa- 
chon (see Childress et al. 201 1, and M. Childress et al. 2012, 
in preparation). The redshifts (from spectroscopic template 
fitting) and morphological types (from visual inspection) are 
listed in Table [T] 

2.3.1. SNIFS Spectrophotometry 

Observations of all six SNe were obtained with SNIFS, 
built and operated by the SNfactory. SNIFS is a fully in- 
tegrated instrument optimized for automated observation of 
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Table 1 

Candidate super-Chandra SN discoveries from SNfactory 



SN Name 


RA 


DEC 


Disc 


UT Date 


Disc. Phase" 


^helio 


Host Type*-' 


E{B-V)mw'' 


SNF 20070528-003 


16:47:31.46 


-1-21:28:33.4 


2007 


May 05 


-7 


0.1171 


dlrr 


0.045 


SNF 20070803-005 


22:26:24.03 


-1-21:14:56.6 


2007 


Aug 03.4 


-11 


0.0315 


Sbc 


0.047 


SN 2007if 


01:10:51.37 


-1-15:27:40.1 


2007 


Aug 25.4 


-10 


0.0742 


dliT 


0.082 


SNF 20070912-000 


00:04:36.76 


-Hl8:09:14.4 


2007 


Sep 12.4 


-17 


0.1231 


Sbc 


0.029 


SNF 20080522-000 


13:36:47.59 


-h05:08:30.4 


2008 


May 22 


-18 


0.0453 


Sb 


0.026 


SNF 20080723-012 


16:16:03.26 


-1-03:03:17.4 


2008 


July 23.4 


-18 


0.0745 


dirr 


0.062 



" In rest-frame days relative to B-band maximum light, as determined from a SALT2 fit to the X'-FS-corrected rest-frame light 
curve. 

From template fit to host galaxy spectrum jChildress et a/.l201 IL M. Childress et a!. 2012, in preparation). 

Moipological type from visual inspection. 

From|Schlegel et al. 1 1 998). 

■= Discovered independently by the Texas Supernova Search <Yuan et a/.ll^DPTlECTTg) and by SNfactory as SNF 20070825-001 
(Scalzo etal. 20 ij). 



point sources on a structured background over the full opti- 
cal window at moderate spectral resolution. It consists of a 
high-throughput wide-ban d pure-lenslet integral field spect ro- 
graph (IFS, "a la TIGER": Bacon et a/.l[T99l l2000[ l200l . a 
multifilter photometric channel to image the field surround- 
ing the IFS for atmospheric transmission monitoring simul- 
taneous with spectroscopy, and an acquisition/guiding chan- 
nel. The IFS possesses a fully filled 6"4 x 6."4 spectroscopic 
field of view (FOV) subdivided into a grid of 15 x 15 spa- 
tial elements (spaxels), a dual-channel spectrograph covering 
3200-5200 A and 5100-10000 A simultaneously, and an in- 
ternal calibration unit (continuum and arc lamps). SNIFS 
is continuously mounted on the south bent Cassegrain port 
of the UH 2.2 m telescope (Mauna Kea) and is operated 
remotely. The SNIFS spectrophotometric da ta reduction 
pipeli ne has been describe d in previous paper s (iBacon et al.\ 
1200 U lAldering et aTll2006t IScalzo ef a/.ll20I0h . We subtract 
the host galaxy ligh t using the methodology described in 
iBongard et al.\ (1201 lb , which uses SNIFS IFU exposures of 
the host taken after each SN has faded away. 

2.3.2. SMARTS Photometry 

The SMARTS imaging of SNF 20070803-005 and 
SN 2007if with ANDICAM consisted of multiple 240 s ex- 
posures in each of the BVRI filters at each epoch. The images 
were proces sed using an automated pip eline based on IRAF 
(lTodyllI993h : this pipeHne was used in IScalzo et aU (120 1 Ol) 
and is described in detail there. We briefly summarize the 
methods below. 

All ANDICAM images were bias-subtracted, overscan- 
subtracted, and flat-fielded by the SMARTS Consortium, also 
using IRAF (ccdproc). In each band, four to six final refer- 
ence images of each SN field were taken at least a year after 
explosion, and combined to form a co-add. This co-add was 
then registered, normalized, and convolved to match the ob- 
serving conditions of each SN image, before subtraction to 
remove the host galaxy light. 

An absolute calibration (zeropoint, extinction and color 
terms) wa s establi s hed o n photometric nights from obser- 
vations of iLandolB (Il992h standards, fitting a zeropoint and 
extinction coefficient for each night separately as well as a 
color term constant across all nights. The calibration was 
transferred to the field stars for each photometric night sep- 
arately using the zeropoint and extinction but ignoring the 
color terms, producing magnitudes on a "natural" ANDICAM 
system which agrees with the Landolt system for stars with 
B-V = V -R = R- 1 = 0. These calibrated magnitudes were 



then averaged over photometric nights to produce final cali- 
brated ANDICAM magnitudes for the field stars. 

Each SN's ANDICAM-system BVRI light curve was then 
measured by comparison to the field stars. To fix the SN's 
location, we found the mean position over observations in 
the same filter, weighting by the signal-to-noise ratio (S/N) 
of each detection. We then measured the final flux in each 
image in a circular aperture centered at this mean location. 
The observer-frame ANDICAM magnitudes were corrected 
for Ga lactic extinction using E( B - V) from ISchlegel et'al] 
(119981) and the extinction law of ICardelU et al.\ (Il988h with 
Ry = 3 . 1 . The SN magnitudes were /T-corrected (iNugent et al.\ 
12002b to rest-frame Bessefl BVRI bandpasses (lBesselllll990l) 
using the SNIFS spectrophotometric time series and a set 
of ANDICAM system throughput curves, the central wave- 
lengths of which were shifted to match ANDICAM observa- 
tions of spectroph otometric standard star s to published syn- 
thetic photometry (IStritzinger et a/.ll2005h . 

2.4. Lightcurves 

The rest-frame. Milky Way de-reddened Bessell BVRI light 
curves of the SNfactory super-Chandra candidates are given 
in Table |2] and shown in Figure [T| The color evolution is 
shown in Figure |2] The light curves of the two SNe at the 
high end of the SNfactory redshift range, SNF 20070528-003 
(z = 0.117) and SNF 20070912-000 (z = 0.123), have less ex- 
tensive coverage than the other SNe; the S/N is lower, and a 
significant fraction of the rest-frame /-band transmission lies 
outside of the observer-frame wavelength range of the SNIFS 
spectrograph. For a small number of observations, rest-frame 
B-band or V-band measurements are unavailable due to in- 
strument problems with the SNIFS blue chan nel. The detai led 
analysis of these lightcurves is presented in ^3.1l and ^3.21 

2.5. Spectra 

Figure [3] presents the subset of spectra that were taken 
near maximum light for the six SNfactory SNe la, including 
SN 2007if. Before maximum light, the spectra show weak 
Si II, S II, and Ca II absorption plus strong Fe II and Fe III ab- 
sorption, in accordance with our selection criteria. SN 199 IT 
and the prototype candidate super-Chandra SN la 2003fg are 
shown for comparison. After maximum light, some features 
noted in SN 2003fg and SN 2007if appear, including the 
iron-peak blend near 4500 A, the s harp n otch near 4130 A 
identified as Cr II in | Scalzo et'oH (120 1 Oh (or tentatively as 
C II A4237 in lHowell et aZ.ll2006l) and fl ie blended Unes nea r 
3300 A identified as Cr II and Co II in IScalzo etaU (I2OI0I) . 
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Table 2 

Rest-frame BVRI light curves 



MJD" Phase'' BVRI Instrument 



SNF 20070528-003 



54250.6 


-7.5 


19.23 ±0.03 


19.30 ±0.03 


19.29 ±0.02 








SNIPS 


54252.5 


-5.8 


19.06 ±0.04 


19.23 ±0.05 


19.23 ±0.03 








SNIPS 


54253.5 


-4.8 


19.00 ±0.03 


19.15±0.04 


19.16±0.03 








SNIPS 


54255.5 


-3.0 


19.00 ±0.02 


19.06 ±0.03 


19.04 ±0.02 








SNIPS 


54261.5 


2.3 


19.04 ±0.02 


18.94 ±0.02 


18.99 ±0.02 








SNIPS 


54263.4 


4.0 


19.15 ±0.02 


19.00 ±0.02 


19.04 ±0.02 








SNIPS 


54268.4 


8.5 


19.34 ±0.02 


19.07 ±0.02 


19.17±0.02 








SNIPS 


54270.4 


10.2 




19.10±0.02 


19.20 ±0.02 








SNIPS 


54270.4 


10.3 




19.10±0.02 


19.22 ±0.02 








SNIPS 


54273.3 


12.9 


19.86 ±0.07 


19.41 ±0.05 


19.47 ±0.04 








SNIPS 


54273.4 


12.9 


19.90 ±0.07 


19.32 ±0.05 


19.36 ±0.04 








SNIPS 


54275.4 


14.7 


20.02 ±0.06 


19.47 ±0.04 


19.46 ±0.03 








SNIPS 


54275.5 


14.8 


19.98 ±0.05 


19.46 ±0.04 


19.45 ±0.03 








SNIPS 


54278.4 


17.4 




19.62±0.12 


19.32 ±0.07 








SNIPS 


54283.4 


21.9 


20.67 ±0.11 


19.89 ±0.05 


19.69 ±0.03 








SNIPS 


54283.4 


21.9 


20.51 ±0.09 


19.88 ±0.05 


19.65 ±0.03 








SNIPS 


SNF 20070803-005 


54318.5 


-8.8 


16.60 ±0.03 


16.68 ±0.03 


16.62 ±0.03 


16 


73±0 


05 


SNIPS 


54320.5 


-6.8 


16.34 ±0.03 


16.38±0.03 


16.35±0.03 


16 


47±0 


05 


SNIPS 


54323.5 


-3.9 


16.10±0.02 


16.12±0.03 


16.13±0.03 


16 


28±0 


05 


SNIPS 


54325.5 


-2.0 


16.10±0.02 


16.11 ±0.03 


16. 15 ±0.03 


16 


34±0 


05 


SNIPS 


54326.8 


-0.7 


16.09 ±0.01 


16.08 ±0.01 


16. 15 ±0.02 


16 


38 ±0 


02 


SMARTS 


54328.7 


1.2 


16.11 ±0.01 


16.06 ±0.01 


16.12±0.02 


16 


39±0 


02 


SMARTS 


54333.5 


5.8 


16.31 ±0.03 


16.13±0.03 


16.17±0.03 


16 


54±0 


06 


SNIPS 


54333.7 


6.0 


16.29 ±0.01 


16.12±0.02 


16.17±0.03 


16 


54±0 


02 


SMARTS 


54335.5 


7.7 


16.40 ±0.03 


16.16±0.02 


16.23 ±0.02 


16 


61±0 


05 


SNIPS 


54336.7 


8.9 


16.46 ±0.02 


16.18±0.02 


16.29 ±0.02 


16 


73±0 


02 


SMARTS 


54338.5 


10.6 


16.67 ±0.03 


16.32 ±0.03 


16.44 ±0.03 


16 


81 ±0 


05 


SNIPS 


54339.7 


11.8 


16.68 ±0.02 


16.33 ±0.02 


16.47 ±0.03 


16 


94±0 


05 


SMARTS 


54340.4 


12.5 


16.78 ±0.03 


16.42 ±0.03 


16.56 ±0.03 


16 


88 ±0 


06 


SNIPS 


54342.7 


14.7 


16.97 ±0.03 


16.51 ±0.04 


16.70 ±0.04 


17 


13±0 


07 


SMARTS 


54343.5 


15.4 


17.10±0.04 


16.62 ±0.04 


16.73 ±0.04 


16 


93±0 


07 


SNIPS 


54345.5 


17.4 


17.32 ±0.04 


16.72 ±0.03 


16.76 ±0.03 


16 


89±0 


05 


SNIPS 


54348.4 


20.2 


17.63 ±0.05 


16.85 ±0.04 


16.78 ±0.04 


16 


80±0 


06 


SNIPS 


54350.4 


22.2 


17.89 ±0.08 


16.96 ±0.06 


16.82 ±0.06 


16 


82±0 


09 


SNIPS 


54353.4 


25.1 


18.20 ±0.09 


17.16±0.05 


16.93 ±0.04 


16 


84 ±0 


06 


SNIPS 


54355.4 


27.0 


18.28±0.13 


17.20 ±0.08 


16.92 ±0.05 


16 


77±0 


07 


SNIPS 


54360.3 


31.8 


18.59 ±0.28 


17.36±0.14 


17.03 ±0.09 


16 


83 ±0 


12 


SNIPS 


54363.3 


34.7 


18.70±0.14 


17.51 ±0.07 


17.19±0.06 


16 


98±0 


08 


SNIPS 


54363.7 


35.1 


18.64 ±0.03 


17.55 ±0.02 


17.22 ±0.02 


17 


19±0 


02 


SMARTS 


54371.7 


42.8 


18.84±0.07 


17.87 ±0.04 


17.59 ±0.04 


17 


58±0 


04 


SMARTS 


54373.3 


44.4 


19.02±0.12 


18.02 ±0.06 


17.76 ±0.06 


17 


59±0 


08 


SNIPS 


54381.6 


52.4 


19.08 ±0.03 


18.21 ±0.02 


18.00 ±0.02 


18 


12±0 


04 


SMARTS 


SNF 20070912-000 


54358.4 


-4.0 


19.30 ±0.02 


19.32 ±0.02 


19.32 ±0.02 








SNIPS 


54360.4 


-2.2 


19.37 ±0.04 


19.27 ±0.05 


19.24 ±0.03 








SNIPS 


54363.4 


0.4 


19.27 ±0.02 


19.16±0.03 


19.13±0.02 








SNIPS 


54365.4 


2.2 


19.29 ±0.03 


19.24 ±0.04 


19.14±0.03 








SNIPS 


54373.4 


9.3 


19.73 ±0.02 


19.42 ±0.02 


19.50 ±0.02 








SNIPS 


54375.4 


11.1 


19.98 ±0.03 


19.57 ±0.02 


19.66 ±0.02 








SNIPS 


54378.3 


13.7 


20.24 ± 0.06 


19.73 ±0.05 


19.77 ±0.04 








SNIPS 


54380.5 


15.6 


20.56 ±0.08 


19.87 ±0.06 


19.79 ±0.04 








SNIPS 


54383.4 


18.2 


20.77 ±0.08 


19.97 ±0.05 


19.89 ±0.03 








SNIPS 


54385.3 


19.9 


20.88 ±0.08 


19.97 ±0.04 


19.70 ±0.03 








SNIPS 


54390.3 


24.4 


21.22±0.11 


20.32 ±0.06 


20.03 ±0.04 








SNIPS 


54390.4 


24.4 




20.34 ±0.06 


20.06 ±0.03 








SNIPS 



" Observer frame JD- 2400000.5. 
In rest-frame days relative to B-band maximum light. 
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Table 2 

Rest-frame BVRI light curves, continued 



MJD" Phase'' BVRI Instrument 



SNF 20080522-000 



54612.4 


-8, 


.7 


17, 


.51 ±0.03 


17.56 ±0.03 


17.54 ±0.03 


17 


.65 ± 0.05 


SNIPS 


54614.4 


-6, 


.8 


17, 


.31 ±0.02 


17.32 ±0.03 


17.34 ±0.03 


17 


.47 ± 0.04 


SNIPS 


54617.3 


-4, 


.0 


17, 


.10±0.02 


17.13±0.03 


17.18±0.03 


17 


.31 ±0.04 


SNIPS 


54619.3 


-2, 


.1 


17, 


.15 ±0.03 


17.16±0.03 


17.20 ±0.03 


17 


.38 ±0.05 


SNIPS 


54622.3 


0, 


.8 


17 


.11 ±0.03 


17.03 ±0.04 


17.05 ±0.03 


17 


.29 ±0.05 


SNIPS 


54624.3 


2, 


.7 


17 


.16±0.03 


17.05 ±0.03 


17.07 ±0.03 


17 


.42 ±0.05 


SNIPS 


54627.3 


5, 


.6 


17 


.24 ± 0.02 


17.10±0.02 


17.13±0.02 


17 


.51 ±0.05 


SNIPS 


54632.3 


10, 


.4 


17 


.48 ±0.03 


17.24 ±0.03 


17.40 ±0.03 


17 


.84 ±0.07 


SNIPS 


54634.3 


12, 


.3 


17 


.66 ±0.03 


17.38±0.03 


17.56 ±0.03 


17 


.95 ±0.07 


SNIPS 


54637.3 


15, 


.2 


18 


.00 ± 0.05 


17.59 ±0.04 


17.74 ±0.04 


18 


.02 ± 0.07 


SNIPS 


54639.3 


17, 


.0 


18 


.25 ±0.05 


17.70 ±0.03 


17.80 ±0.04 


18 


.02 ±0.07 


SNIPS 


54642.3 


19, 


.9 


18 


.59 ±0.07 


17.83 ±0.04 


17.80 ±0.04 


17 


.94 ± 0.07 


SNIPS 


54644.3 


21 


.8 


18 


.76 ±0.09 


17.90 ±0.05 


17.82 ±0.04 


17 


.89 ±0.08 


SNIPS 


54647.3 


24, 


.7 








17.84 ±0.05 


17 


.80 ±0.08 


SNIPS 


54650.3 


27, 


.6 


19 


.23 ±0.09 


18.13±0.04 


17.88 ±0.03 


17 


.75 ± 0.05 


SNIPS 


54652.3 


29, 


.5 


19 


.46 ±0.14 


18.23 ±0.05 


17.93 ±0.04 


17 


.76 ±0.06 


SNIPS 


54654.3 


31 


.4 


19 


.55 ±0.16 


18.38 ±0.06 


18.07 ±0.05 


17 


.87 ±0.07 


SNIPS 


54662.3 


39, 


.1 


19 


.88 ±0.25 


18.75 ±0.09 


18.46 ±0.07 


18 


.22 ±0.09 


SNIPS 


54668.3 


44, 


.8 


20 


.01 ±0.16 


18.98 ±0.07 


18.74 ±0.07 


18 


.59±0.10 


SNIPS 


54669.3 


45, 


.7 


19 


.92 ±0.43 


19.05±0.17 


18.86±0.14 


18 


.70±0.21 


SNIPS 


54672.3 


48, 


.6 


20, 


.23 ±0.30 


19.10±0.13 


18.89±0.11 


18 


.88±0.17 


SNIPS 


SNF 20080723-012 


54674.3 


-6, 


.2 


18, 


.07 ±0.02 


18.18±0.02 


18.11 ±0.02 


18 


.23 ±0.03 


SNIPS 


54677.3 


-3, 


.4 


18, 


.04 ±0.03 


18.12±0.04 


18.07 ±0.03 


18 


.17 ±0.05 


SNIPS 


54679.3 


-1 


.5 


18, 


.00 ±0.02 


17.96 ±0.02 


17.96 ±0.02 


18 


.16±0.03 


SNIPS 


54682.3 


1 


.2 


17 


.96 ±0.02 


17.92 ±0.03 


17.94 ±0.02 


18 


.19 ±0.04 


SNIPS 


54684.4 


3, 


.1 


18 


.08 ±0.06 


17.97 ±0.06 


17.85 ±0.05 


17 


.67 ± 0.09 


SNIPS 


54687.3 


5, 


.9 


18 


.12 ±0.02 


17.88 ±0.02 


17.93 ±0.02 


18 


.26 ± 0.04 


SNIPS 


54689.3 


7, 


.7 


18 


.31 ±0.02 


17.99 ±0.02 


18.07 ±0.02 


18 


.46 ± 0.04 


SNIPS 


54694.3 


12, 


.4 


18 


.92 ±0.10 


18.38±0.07 


18.44 ±0.05 


18 


.69 ±0.10 


SNIPS 


54696.3 


14, 


.2 


19 


.03 ±0.08 


18.42 ±0.05 


18.42 ±0.03 


18 


.58 ±0.07 


SNIPS 


54699.3 


17, 


.0 


19 


.36 ± 0.05 


18.56±0.03 


18.47 ±0.03 


18 


.52 ± 0.05 


SNIPS 


54702.3 


19, 


.8 


19 


.65 ± 0.04 


18.71 ±0.03 


18.55 ±0.03 


18 


.54 ± 0.05 


SNIPS 


54704.3 


21 


.7 


19 


.97±0.18 


19.04 ±0.09 


18.78 ±0.06 


18 


.47 ±0.10 


SNIPS 


54707.3 


24, 


.5 


20 


.15±0.14 


19.01 ±0.06 


18.75 ±0.04 


18 


.57 ±0.07 


SNIPS 


54709.3 


26, 


.3 


19 


.97 ±0.09 


18.96 ±0.04 


18.69 ±0.03 


18 


.58 ±0.05 


SNIPS 


54712.3 


29, 


.1 


20 


.16 ±0.09 


19.11 ±0.04 


18.82 ±0.04 


18 


.67 ±0.06 


SNIPS 


54714.3 


31 


.0 


20 


.34 ±0.10 


19.23 ±0.05 


18.93 ±0.04 


18 


.76 ±0.06 


SNIPS 


54719.3 


35, 


.6 


20 


.39±0.12 


19.42 ±0.06 


19.13±0.05 


18 


.93 ± 0.08 


SNIPS 


54726.3 


42, 


.1 


20 


.67 ±0.15 


19.66 ±0.09 


19.39 ±0.08 


19 


.18±0.12 


SNIPS 


54734.2 


49, 


.6 


20 


.46 ±0.22 


19.87±0.16 


19.71 ±0.13 


19 


.61 ±0.22 


SNIPS 



" Observer frame JD- 2400000.5. 
In rest-frame days relative to B-band maximum light. 
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Figure 1. Rest-frame Bessell BVRI light curves from ANDICAM+SNIFS 
from SNIFS flux-calibrated spectroscopy. Circles: ANDICAM BVRI, K 
red lines mark the light curve phase range corresponding to the break in v 

(The 4130 A feature is not clearly present in the z = 0.123 
candidate SNF 20070912-000, although thi s may be due to 
the lower S/N of the spectrum.) Taubenbe rger et al.\ (1201 lb 
note that a similar feature in SN 2009dc strengthens with time 
rather than fading, as one would expect for Cr II rather than 
Cll. 

SN 2007if also shows a weak C II A6580 li ne in the post- 
maximum spectra, which IScalzo et al.\ (120101) interpreted as 
a signature of unburned material from the explosion. While 



-10 10 20 30 40 50 60 
Rest-Frame Days Since B-Band Maximum 

Inverted triangles: Rest-frame Bessell BVRI magnitudes synthesized 
corrected to the re spec tive rest-frame Bessell filters. Vertical dotted 
(Si II)discussed in 33.31 



this line is not detected unambiguously in the other SNe, the 
unusually shallow slope of the red wing of Si II A6355 e.g. 
in SNF 20070528-003 and SNF 2( )080723-012 inay be a sign 
that C II is present in these SNe (iThomas ef a/.ll20il1). T he 
spectral properties of the full dataset are analyzed in 



3. ANALYSIS 



The analysis (this section), modeling (® and interpreta- 
tion (^ for our sa mple parallels that made for SN 2007if in 
IScalzo et al.\ (120101) . with improvements described below. We 



Scalzo et al. 



SNF 20070528-003 




SNF 20070803-005 




SN 2007if 



SNF 20070912-000 



3 

2 

o 

o 1 
U 



-1 



T 

V 



V 
T 




SNF 20080522-000 



SNF 20080723-012 



3h 
2 



O 

o 1 
U 



- 



T" 



T" 



Vv VV W V 



v7 



-10 10 20 30 40 50 60 
Rest-Frame Days Since B-band Maximum 




-10 10 20 30 40 50 60 
Rest-Frame Days Since B-band Maximum 



Figure 2. Rest-frame Bessell color evolution from ANDICAM (circles) and SNIFS spectrophotometry (inverted triangles). Dotted lines indicate 
fits to the Lira relation with a floating excess. 



include a re-analysis of our observations of SN 2007if in this 
paper using the improved techniques, for a more direct com- 
parison with the other SNe presented in this paper. 

3.1. Maximum-Light Behavior, Colors and Extinction 

As discussed in IScalzo et a/1 (I2010I) . SN 2007if shows 
no distinct second maximum. While Figure [T] does show 
second maxima for our other SNe, its prominence is sup- 
pressed relative to normal SNe la. In SNF 20070803-005 and 
SNF 20080723-012, the peak-to-trough difference in a quin- 



tic polynomial fit to the data from day +7 to day +42 is only 
0. 1 1 mag, vs. 0.23 mag for SNF 200805 22-000 a nd 0.3 5 mag 
for the SALT2 model with a1 = 1, c = 0. iKasenI (120061) noted 
three physical effects which could reduce the contrast of the 
/-band second maximum; low ^''Ni mass, efficient mixing of 
""^Ni into the outer layers of ejecta, and greater absorption in 
the Ca II NIR triplet line source function. Since our SNe are 
all overluminous with broad light curves, Arnett's rule gives 
a high nickel mass, as we find in the next section. The /-band 
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Figure 3. Selected sp ectra near m aximum light of SNfac tory candidate super-Chandra SNe la (black), with spe ctra of SN 199 IT ( red) at -6 d 
dMazzali ef"allll995l) and at +7 d jFilippenko et fl/.lll992l) . and the candidate super-Chandra events SN 2003fg dHowell er alll200d green) and 
SN 2007if dScalzo et a/.llMol , green) for comparison. Spectra have been smoothed with a 15-point (width 2000 km s ') Savitzsky-Golay filter 
for presentation purposes. Consecutive post-maximum spectra of the z ~ 0.12 S Ne SNF 20070528-0 03 and SNF 20070912-000 have been 
co-added to improve signal-to-noise. Common features identified with SYNAPPS dThomas et fl/.l20l"l1) are marked. 
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Table 3 

Derived quantities from SALT2 fits to light curves 



SNName MJD(B,„ax) Mv" .vi c y^jv l\iini(Bf (B-V)n,ax Bmax-Vmax i^mif 



SNF 20070528-003 


54258.7 


-19.71 ±0.02 


1.08 ±0.20 


0.02 ±0.03 


1.82 


0.84 ±0.05 


0.01 ±0.04 


0.03 ±0.04 


-0.48 


SNF 20070803-005 


54327.6 


-19.59 ±0.04 


1.04 ±0.09 


0.01 ±0.03 


0.74 


0.85 ±0.05 


0.00 ±0.04 


0.02 ±0.04 


-0.28 


SN 2007if 


54346.8 


-20.38 ±0.02 


1.41±0.13 


0.12±0.03 


3.48 


0.79 ±0.06 


0.12±0.05 


0.14 ±0.04 


-1.34 


SNF 20070912-000 


54363.0 


-19.69 ±0.02 


-0.03 ±0.19 


0.04 ±0.03 


1.41 


1.08 ±0.04 


0.05 ±0.02 


0.06 ±0.02 


-0.52 


SNF 20080522-000 


54621.4 


-19.59 ±0.03 


1.18±0.13 


-0.01 ±0.03 


0.43 


0.83 ±0.06 


-0.02 ±0.04 


0.00 ±0.04 


-0.17 


SNF 20080723-012 


54680.9 


-19.71 ±0.02 


0.53±0.12 


0.04 ±0.03 


2.61 


0.93 ±0.04 


0.03 ±0.03 


0.05 ±0.03 


-0.59 



" Evaluated from SALT2 rest-frame Vmax, with distance modulus at the host galaxy redshift assuming a ACDM cosmology with SIm = 0.28, fiA = 0.72, 
//o = 72 km s"' Mpc"' . 

^ Light curve decline rate, evaluated directly from the best-fit SALT2 mod el, accounting for erro r in the date of B-band maximum light. 
Hubble residuals from ACDM cosmology, evaluated using equation 2 of lSuUivan et ali 1201 II) . 
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first maximum is roughly concurrent with B-band maximum 
for our S Ne, rather tha n being significantly delayed (see fig- 
ure 14 of lKasenll2006l) . so it seems unlikely that emission in 
the Ca II NIR triplet is contributing significantly. The most 
likely interpretation, especially given the prominence of Fe- 
peak lines in early spectra of our SNe, is that ^*Ni is well- 
mixed into the outer layers^ 

Following practice from lScalzo et a/.l (120101) . we use the up- 
dated version (v2.2) of the SALT2 light curve fitter jGuv et al\ 
120101) to interpolate the magnitudes and colors of each SN 
around maximum light, and to establish a date of Z?-band 
maximum with respect to which we can measure light-curve 
phase. While we use SALT2 here as a convenient functional 
form for describing the shape of the light curves near maxi- 
mum light, and to extract the usual parameters describing the 
light curve shape, we do not expect the SALT2 model, trained 
on normal SNe la, to give robust predictions for these peculiar 
SNe la outside the phase and wavelength coverage for each 
SN. To minimize the impact of details of the SALT2 spectral 
model on the outcome, we use SALT2 in the rest frame, in- 
clude SALT2 light curve model errors in the fitting, and we fit 
BVR bands only; /-band is excluded from the fit. The quanti- 
ties derived from the SALT2 light curve fits are shown in Ta- 
blets] A cross-check in which cubic polynomials were fitted to 
each band produces peak magnitudes and dates of maximum 
in each band consistent with the SALT2 answers, within the 
errors, for all of the new SNe; we adopt the SALT2 values as 
our fiducial values for direct comparison with other work. 

We estimate the host reddening of the SNe in two sepa- 
rate ways. First, we fit the B-V color behavior of each 
SN to the Lira relation (iPhillips et al\ 119991: iFolatelli etal\ 
l2009h for those SNe for which we have appropriate light 
curve phase coverage. The Lira relation is believed explicitly 
not to hold for SN 2007if and the candidate super-Chandra 
SN l a 2009dc ( Yamanaka ef a/. 2009; Taubenberg er et al\ 
1201 ll) . but its value may nevertheless be useful in study- 
ing the relative intrinsic color of these SNe. Additionally, 
we search for Na I D absorption at the redshift of the host 
galaxy for each SN. We perform a X' fit to the Na I D 
line profile, modeled as two separate Gaussian lines with 
full width at half maximum equal to the SNIFS instrumen- 
tal resolution of 6 A, to all SNIPS spectra of each SN, as for 
SN 2007if in IScalzo et all (120101) . In the fit, the equivalent 
width EWCNa I D) of the Na I D line is constrained to be non- 
negative. We convert these to estimates of £'(B-V')host using 
both the shall ow-slope (0.16 A~') and steep-slope (0.51 A"') 
relations from lTuratto. Benetti & Cappellard (12002b ("TBC"). 
While the precision of these relatio ns has been called int o 
question when used on their own (e.g. lPoznanski et aZll201 Ih . 
we believe that examining such estimates together with the 
Lira relation and fitted colors from the light curve can pro- 
vide helpful constraints on the importance of host reddening. 
The best-fit Lira excesses, values of ^^(Na I D), and derived 
constraints on the host galaxy reddening are listed in Table ID 

We detect weak Na I D absorption in SNF 20070803-005 
and SNF 20080522-000. Neither of these SNe appear to 
have very red colors according to the SALT2 fits, and we 
believe it to be unlikely that either are heavily extinguished, 
so for purposes of extinction corrections to occur later in our 
analysis, we use reddening estimates from the shallow-slope 
TBC relation, together with a CCM dust law with Ry = 3.1 
dCardelli et al. 19883- When appUed to Milky Way Na I D 
absorption in our spectra, the shallow-slope TBC relation 



Table 4 

Host galaxy reddening estimates from the Lira relation and £'W(Na I D) 

fits 





NalD 


Shallow 


Steep 


Lira 


SN Name 


EW" (A) 


TBC'' 


TBC^ 


Relation'' 


SNF 20()70528-()03 
SNF 20070803-005 
SN 2007if 
SNF 20070912-000 


00+" 
00+0- '2 

"■""-0.00 
"■^^-0.18 


00+0 05 

"■""-0.00 

02+0 01 

"■"^-0.0 1 

00+0 01 

"■""-0.00 

04+0 07 

"■"^-0.03 


00+0- '2 
"■""-0.00 

06+^ 03 

"■""-0.03 

00+* 02 

"■""-0.00 

12+*^22 
"■^ -0.09 


0.06 ±0.06 
0.13 ±0.07 



— u. 1 a —\i.y)D — u-Ui" 

SNF 20080522-000 0.29!j;;[]^ ^■^^'^am O.ll^jS g^ 0.13 ±0.12 
SNF 20080723-012 OW*^^^ "■"O^mo "■"O^oioo 0.00±0.10 



" Measured from a simultaneous fit of the Na I D absorption line profile 
to all SNIFS spectra of each SN. 

E(B — 'V) derived from Na I D absorption, using the "shallow" slope 
(0.16 A~') of lTuratto. Benetti & Cappellard <2D02|) . 

E{B - V) derived from Na I D absorption, using the "steep" slope 
(0.51 A"') of lTuratto. Benetti & Cappellard <2D02|) . 
** Best-fit E(B — V ) from the Lira relation in the form given in 
jPhillips e? a/.lll999l) . EiTor bars include a 0.06 mag intrinsic dispersion 
of normal SNe la around the relation, added in quadrature to the statisti- 
cal errors. 

P roduc es E{B-V) estimates consistent with ISchlegel et dl\ 
19981) . SNF 20070912-000 shows a marginal (< 2a) detec- 
tion, though the spectra are noisy and the limits are not strong. 
We detect no Na I D absorption in the other SNe. 

Based on the very strong limit on Na I D absorption from 
the host galaxy, IScalzo et al\ ( 120101) inferred that the large 
Lira excess of SN 2007if was not due to host galaxy extinc- 
tion. Those SNe observed at sufficiently late phases show 
measured Lira excesses consistent with zero, additional ev- 
idence that host galaxy dust extinction is minimal for these 
SNe if the Lira relation holds. 

Allowing for varying amounts of extinction associated with 
the Lira excess or Na I D absorption, each of the new SNe 
have maximum-light {B - V) colors consistent with zero. 
While most of our SNe have well-sampled light curves around 
maximum light and hence have well-measured maximum- 
light colors, SN 2007if has a gap between -9 days and 
+5 days with respect to B-band max imum (phases fixe d by 
the SALT2 fit). The light curve fit of (IScalzo ef a/.ll2010l) . us- 
ing SALT2 v2.0 and a spectrophotometric reduction using a 
previous version of the SNIFS pipeline, suggests a red color 
B-V = 0. 16 ±0.06. The more recent reduction fit with SALT2 
v2.2 gives B-V = 0.12 ±0.06. However, the directly mea- 
sured B-V color of SN 2007if at -9 days (-0.07 ± 0.04) and 
at +5 days (0.15 ±0.02) are each consistent within the er- 
rors with the mean values at those epochs for our other five 
SNe. The systematic error on the maximum-light color of 
SN 2007if, at least 0.04 mag, may therefore be too large for 
it to be considered significantly redder at maximum than its 
counterparts. 

3.2. Bolometric Light Curve and ^^Ni Synthesis 

As input to our further analysis to calculate ^''Ni masses 
and total ejected masses for our sample of SNe, we calcu- 
late quasi-bolometric UVOIR light curves from the photom- 
etry. To derive bolometric fluxes from SNIFS spectropho- 
tometry, we first der edden the spectra to account for Milky 
Way dust reddening (ISchlegel et a/.|[T998l) . then deredden by 
an additional factor corresponding to a possible value of the 
host galaxy reddening, creating a suite of spectra covering the 
range 0.00 < E{B-V\o^i < 0.40 in 0.01 mag increments. We 
then integrate the dereddened, deredshifted, flux-calibrated 
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spectra over all rest-frame wavelengths from 3200-9000 A. 
For each quasi-simultaneous set of ANDICAM BVRI obser- 
vations and each possible value of E{B -V)host, we multiply 
the dereddened, deredshifted SNIFS spectrum nearest in time 
by a cubic polynomial, fitted so that the synthetic photom- 
etry from the resulting spectrum matches, in a least-squares 
sense, the ANDICAM imaging photometry in each band. We 
then integrate this "warped" spectrum to produce the bolo- 
metric flux from ANDICAM. This procedure creates a set of 
bolometric light curves with different host galaxy reddening 
values which we can use in our later analysis (see ® . 

SNF 20070528-003 and SNF 20070912-000 are at a higher 
redshift than our other SNe, such that SNIFS covers only the 
rest-frame wavelength range 3000-8500 A, so we integrate 
their spectra in this range instead. We expect minimal sys- 
tematic error from the mismatch, since the phase coverage for 
these two SNe is such that only the bolometric flux near max- 
imum light, when the SNe are still relatively blue, is useful 
for the modeling described in ^ 

To account for the reprocessing of optical flux into the near- 
infrared (NIR) by iron-peak elements over the evolution of 
the light curve, we must apply a NIR correction to the in- 
tegrated SNIFS fluxes. No NIR data were taken for any of 
the SNe except SN 2007if; in general they were too faint 
to observe effectively with the NIR channel of ANDICAM. 
Since these are peculiar SNe, any correction for the NIR flux 
necessarily involves an extrapolation. Since the /-band sec- 
ond maximum has low contrast for all the SNe in our sam- 
ple, the JHK behavior should be similar for th ese SNe to th e 
extent that / and JHK are related (e.g. as in iKasenI 120061) . 
With these caveats, we therefore use the NIR corrections for 
SN 2007if derived in lScalzo ef all JMTol) for all of our SNe, 
with the time axis stretche d according to th e stretch factor de- 
rived from the SALT2 xi (iGuy et al.\ 2007h to account for the 
different timescales for the development of line blanketing in 
these SNe. Our modeling requires knowledge of the bolo- 
metric flux only near maximum light (to constrain the ^^Ni 
mass) and more than 40 days after maximum light (to con- 
strain the total ejected mass). The NIR correction is at a min- 
imum near maximum light (~ 5%), and at a maximum near 
phase +40 days 25%), so it should not evolve quickly at 
these times and our results should not be strongly affected. 
We assign a systematic error of ±5% of the total bolometric 
flux (or about ±30% of the NIR flux itself near +40 d) to this 
correction. 

To estimate ^^Ni masses, we also need a measurement of 
the bolometric rise time frise.boi- We establish the time of bolo- 
metric maximum light fmax.boi by fitting a cubic polynomial to 
the bolometric fluxes in the phase range -10 d < f < +20 d. 
We then calculate the final bolometric rise time frise.boi via 



Table 5 

^^Ni mass reconstruction for the SNfactory SNe 



frise.boi — frise.S fmax,S "t" fmax.bol- 



(1) 



We find fmax.boi to occur about 1 day earlier than tmax.B for the 
SNe in our sample. We have a strong constraint on the time of 
explosion only for SN 2007if (IScalzo ef a/.ll2010h . for which 
this procedure results in a bolometric rise time of 23 days. We 
have utilized the discovery data from our search along with 
our spectrophotometry to constrain the rise times of the new 
SNe presented here. For these we find f„.se,B = 21 .2 ± 1 .9 days, 
and so use trise.hoi = 20 ± 2 days in our models. Our value is 
very similar to the value of frise b = 21 zfc 2 days given by th e 
sample of 1991T-like SNe la in lOaneshaUngam et al.\ ( 1201 Ih . 
This approach leads to more conservative uncertainties than 



SN Name 


Lbol (lO'*^ erg 




'risc.bol (d)"" 


Mn, 


(Mq)'' 


SNF 20070528-003 


1.87 ± 


0.06 


20 ±2 


0.78 


±0.11 


SNF 20070803-005 


1.62 ± 


0.05 


20±2 


0.69 


±0.07 


SN 2007if 


2.93 ± 


0.22 


23 ±2 


1.38 


±0.09 


SNF 20070912-000 


1.86 ± 


0.06 


20±2 


0.77 


±0.10 


SNF 20080522-000 


1.70 ± 


0.06 


20±2 


0.74 


±0.08 


SNF 20080723-012 


1.82 ± 


0.07 


20±2 


0.76 


±0.10 



" Calculated from Equation[T] 
Assuming fiducial a = 1.3. 



the single-stretc h correc t ion oflConlev et al.\ (120061) : Figure 6 
of Ganeshalinga m ef all ( 1201 Ih suggests that the relation be- 
tween rise time and decline rate Amis (or stretch s) may break 
down at the high-stretch end. 

We calculate the ^^Ni mass, M^i, for the six SNe by relating 
the maximum-light bolometric l uminosity Lboi to the luminos- 
ity from radioactive decay Ljad (lArnettlll982) : 



ad 



56NiA56Ni 



1,7 

-^56Co 



A56Ni - A56C0 



X (g-^'OCo'-g-^^ONi') , (2) 

where t (= fi-ise.boi in this case) is the time since explosion, 
A?56Ni = M56Ni/(56 AMU) is the number of ^^Ni atoms pro- 
duced in the explosion, Ass^i and Xsico ^re the decay con- 
stants for ^^Ni and ^^Co (e-folding lifetimes 8.8 days and 
111.1 days) respectively, and QssNi^, Qaco.-/ 25'iCo.e+ 
are the en ergies released in the different stages of the de- 
cay chain (iNadvozhinI 1 1 9941) . The dimensionless number a 
is a correction factor accounting for the diffusion time de- 
lay of gamma-ray energy through the ejecta, typically rang- 
ing between 0.8 and 1.6 for reasonable explosion models (see 
e.g.. Table 2 of Jioflich & Khohkloy 1996|, where it is called 
0. A nom inal value of a = 1 . 2 is often used in the lit- 
erature (e.g. iNugent et al.\ 1 1 9951: [Branch & Khokhlovl [T99l 
Howell et a/.ll2006r 20091) . The tamped detonation mode ls of 
iKhokhlov ef g/.l ( Il993h and iHoflich & Khohklovl ([T996h . on 
which we will base our modeling later in the paper, have 
slightly higher values closer to a = 1.3. We therefore adopt 
a fiducial value of a = 1.3 for our simple estimate here. 

The resulting ^^Ni mass estimates are shown in Table|5] The 
new SNe have MNi in the range 0.7-0.8 M©, at the high end 
of what might be expected for Chandrasekh ar-mass explo- 
sions: the well-known W7 deflagration model (iNomoto et al.\ 
[19841) produced 0.6 Mq of ^*Ni, while some delayed detona- 
ti on models can produce up t o 0.8 M© (e.g., the N21 model 
of lHoflich & Khohklovll996l) . 

3.3. Spectral Features and Velocity Evolution 

The time evolution of the position of the Si II A6355 ab- 
sorption minimum, showing the recession of the photosphere 
through the ejecta, is shown in Figure|4l The measurements of 
the absorption minimum were made as follows: bins in each 
spectrum immediately to the right and left of the line feature 
were used to fit a linear pseudocontinuum, ^A.cont — ^ 

+ bX. 

This fitted pseudocontinuum was divided into each bin of the 
spectrum in the region of the line. The resulting spectrum was 
smoothed with a third-order Savitzsky-Golay filter and the 
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Figure 4. Time evolution of the velocities of various intermediate-mass-element absorption minima from the super-Chandra candidate sample: 
Si II A6355, S II A5640, and Ca II H-l-K. The thick dash-dot line shows a fit to a constant for all i '(Si II) data before th e end of the plateau 
phase for each SN. The thin dashed lines show the 1-sigma range of behavior for the LVG subclass of iBenetti et al.\ ( 120051) . 
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Table 6 

Velocity gradient cliaracteristics 



5800 6000 6200 6400 
Rest-Frame Wavelength (A) 



6600 



Figure 5. Detailed evolution of SNF 20070803-005 spectra around 
the Si II A6355 feature. 

minimum was recorded as the bin with the lowest signal. The 
error bars on the procedure were determined through a boot- 
strap Monte Carlo: In the first stage, values of a and b rep- 
resenting possible pseudocontinua /^x,cont were sampled using 
the covariance matrix of the pseudocontinuum fit; for each 
candidate pseudocontinuum, fluctuations typical of the mea- 
sured errors on each spectral bin were added to the smoothed 
spectrum, and the results were smoothed and the minimum 
measured again. We have verified that the Savitzsky-Golay 
filter preserves the line minimum, so that smoothing a spec- 
trum twice introduces negligible systematic error. The final 
velocity values and their errors were measured as the mean 
and standard deviation of the distribution of absorption mini- 
mum velocities thus generated. While this method has slightly 
less statistical powe r than a fit to the entire line profile as in 
IScalzo et al\ (l20Tol) . we believe it is more robust to possible 
bias from line profiles with unusual shapes, and provides more 
realistic error bars for the line minimum. 

We found that for lines with equivalent widths less than 
15 A, the absorption minima had unreasonably large uncer- 
tainties and/or showed large systematic deviations from the 
trend described by stronger measurements. When measuring 
these absorption minima, we are proba bly simply me a suring 
uncertainty in the pseudocontinuum. iBlondin et al.\ (1201 Ih 
saw similar effects when measuring velocities of very weak 
absorption minima, to the extent that the Si II A6355 velocity 
would even be seen to increase with time. We therefore reject 
measurements of such weak absorption features. 

Our candidate super-Chandra SNe la share a slow evolu- 
tion of the Si II velocity, consistent within the errors with 
being constant in time for each SN from the earliest phases 
for which measurements are available. Table |6] shows the 
fitted constant velocities and the chi-square per degree of 
freedom, xt^ for a fit to a constant. For comparison with 
earlier work, the velocity gradient v = -dv/dt calculated as 



SN Name 






1 


;b 




At/ 




SNF 20070528-003 


9371 


± 


171 


-6 


± 


29 


-4 


16 


0.85 


SNF 20070803-005 


9695 


± 


81 


46 


± 


23 


-1 


14 


1.00 


SN 2007if 


8963 


± 


248 


31 


± 


29 


-9 


21 


0.76 


SNF 20070912-000 


9201 


± 


403 














SNF 20080522-000 


10936 




107 


5 




10 


-4 


16 


0.41 


SNF 20080723-012 


10391 




291 


-72 




47 


-2 


10 


0.72 



" Best-fit constant velocity, in km s ' , with eiTor. 

^ Slope of the best-fit line to absorption line velocities from the first re- 
liable measurement until the break associated with Fe II line blending, in 
km s"' day"' . 

Phase of first available measurement, in days with respect to B-band 
maximum light, which we interpret to be the start of the plateau. 

Minimum duration of the visible plateau phase in days, until Si II be- 
comes blended with Fe II. 

Chi-square per degree of freedom for fit to a constant. 



the slope of the best-fit linear trend of the measurements 
before day -1-14, is also listed, along with the formal error 
from the fit . All of our SNe would be classified as Benetti 
LVG ( Benetti et a/.l 12005 ) based on their velocity gradients. 
While the slope of the straight-line fit to the absorption ve- 
locities for SNF 20070803-005 seems to differ from zero at 
the 2(7 level (formal errors), the reduced chi-square for this 
fit is extremely small {xl = 0.006) and we conclude that the 
evolution cannot be reliably distinguished from a constant 
ixi = 1.00). Similarly, the best-fit line to the absorption ve- 
locities for SNF 20080723-012 has a positive slope at 1. 5a 
(v = 72±47 km s~' day"'), but once again, a constant is a 
good fit to the data (xl = 0.72). 

The slow Si II absorption velocity evolution usually lasts 
until about two weeks after maximum light, after which a 
break in the behavior occurs and a pronounced decline begins, 
at a rate of ^ 500 km s~' day"' . At this point, developing Fe II 
lines have probably blended with S i II and made the ve locity 
measurements unreliable (see, e.g. iPhillips et a/.lll992h . The 
transition occurs concurrently with the onset of the second 
maximum in the /-band light curve, which may be attributed 
to light reprocessed from bluer wavelengths by the re combi- 
nation of Fe III to Fe II as the ejecta expand and cool dKasenl 
120061) . Figure |5] shows the evolution of the Si II A6355 line 
profile for SNF 20070803-005 near the transition. The line 
profile after the break shows a portion of the Si II line near the 
plateau velocity, suggesting that the material responsible for 
the plateau has thinned but is not yet transparent. We mark 
the phases of the last v(Si II) measurement visually consistent 
with the plateau velocity, and the first measurement inconsis- 
tent with it, by vertical dashed lines in Figure [T] to show their 
correspondence with the /-band second maximum. 

In four of the six SNe, the plateau velocity is low 
(~ 9000 km s"'), inconsistent with the normal ran ge of be- 
havior of the LVG subclass of iBenetti et al.\ (120051) . The re- 
maining two SNe, SNF 20080522-000 and SNF 20080723- 
012, have higher plateau velocities (~ 1 1000 km s"'), falling 
roughly into the range of LVG behavior, although their veloc- 
ity gradients are still flatter than any remarked upon in that 
work. Due to the lower S/N of our spectra of SNF 20070912- 
000 and the weakness of the Si II A6355 absorption feature, 
only two measurements of the line velocity before the break, 
each around 9000 km/s with large errors, were extracted; we 
can at most say that the observed velocity gradient in this SN 
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is consistent with that of the other SNe in our sample. 

Figure |4] also shows the line minimum velocities of 
Ca II H+K and S II AA5456,5640. Although, as expected, 
Ca II H+K stays optically thick longer than Si II A6355, its ve- 
locities at early times are consistent with the observed plateau 
behavior shown by Si II. The S II lines are weak, but when 
their absorption minima can be reliably measured, they do not 
show dramatically higher or lower velocities than the other 
lines. This supports the interpretation that all of these lines 
are formed in the same thin, dense layer of ejecta. 

4. CONSTRAINTS ON TOTAL MASS AND DENSITY STRUCTURE 

The modeling proc edure we use here represents a refined 
version of that used in lScalzo et g/.l (120101) . w hich we compare 
and contrast with the similar approach of iStritzinger et al.\ 
(|2006) in gJjfor the special case of a simple equivalent ex- 
ponential density profile. We then describe our extensions to 
t he m ethod, including modeling of density profiles with shells 
04.2l i. priors on the central density 04.3l l, and calculation of 
the ^^Ni form factor q 04.4l i. We then present our final mod- 
eling results in ^4.51 

4.1. SN la Ejected Mass Measurements Using the Equivalent 
Exponential Density Formalism 

The ejecta density structure of SNe la is frequently modeled 
as an exponential p(v) cx exp(-v/ve), where the ejecta are in 
homologous expansion at velocity v = (r- rinitiaO/f since the 
explosion at time f = 0, and is a characteristic velocity scale. 
Many hydrodynamic mod els of SN la explosio ns, including 
the well-known W7 model ( iNomoto et a/.ll984h . have density 
p rofiles which ar e very close to exponential. 

iJefferyl (Il999h made semianalytic calculations of the time 
evolution of the gamma-ray energy deposition in an exponen- 
tial model SN la, with reference to its bolometric light curve. 
By about 60 days after explosion, virtually all of the ^^Ni has 
decayed and the dominant energy source is the decay of ^^Co 
(which in turn was produced by ^''Ni decay at earlier times). 
The optical depth to Compton scattering of ^^Co gamma rays 
behaves as = (f/fo)~^ with = 1 at some fiducial time fo- 
The value of fo can be extracted by fitting the bolometric light 
curve for f > 60 d to a modified version of Equation |2] in 
which a = 1 (rather than its maximum-light value from Ar- 
nett's rule) and the term corresponding to ""^Co gamma rays 
is multiplied by a factor 1 -exp(-T-y). The total mass of the 
ejecta can then be expressed as 

MwD=— (vA)^ (3) 

Here, is the Compton scattering opacity for ^^Co gamma 
rays, and ^ is a form factor describing the distribution of the 
^^Co in the ejecta, which follows the original distribution of 
^^Ni in the explosion. The value of is expected to lie in 

the range 0.025-0.033 cm"- g (ISwartz ef a/.ll 19951) . with the 
low end (0.025) corresponding to the optically thin regime. 
The value of q can be readily calculated given an assumed 
distribution of ^''Ni (see gjUbelow). 

Stritzing er et al\ (|2006) used this method to measure pro- 
genitor masses for a sample of well-observed SNe la with 
UBVRI light curve coverage. They constructed "quasi- 
bolometric" Ught curves according to the procedure of 
IContardo. Leibundgut & Vaccal ( 120001) . by converting the ob- 
served UBVRI magnitudes to monochromatic fluxes at the 



central wavelengths of their respective filters, then summing 
them, using corrections for lost flux between filters derived 
from spectroscopy of SN 1992 A. Th ey then applied the semi- 
analytic approach of iJeffervl ( 119991) to fit for the gamma-ray 
escape fraction, and hence the ejected mass of the SN la pro- 
genitor. 

Our own work improves on previous use of this method in 
two important ways. First, IStritzinger et al\ (120061) made no 
attempt to correct for the NIR contribution to the bolometric 
flux, simply asserting that it is small during the epochs of in- 
terest. In Scalzo et al.\ (1201 Oh we found that for SN 2007if the 
NIR contribution was indeed small (~ 5%) near maximum 
light, but was greater than 25% at 40 days after maximum 
light, and our estimate for its value at 100 days after maxi- 
mum light is still around 10%. Therefore, at least for SNe la 
like the ones we study here, the method of IStritzinger et ai\ 
( 120061) underestimates the fraction of trapped ^*Co gamma- 
rays for a given initial ^^Ni mass, and hence the ejected mass, 
as a result of neglecting NIR flux. 

Second, our fitting procedure includes covariances between 
different inputs to the prediction for the bolometric light 
curve, constrained by a set of Bayesian priors motivated by 
explosion physics. Specifically, covariances between q, Vg, 
and a may influence the interpretation of the fitted value of fo. 
IStritzinger et al\ (120061) simply fixed the ^''Ni mass from Ar- 
nett's rule, and then fit for fo. Similarly, they assume q = 1/3, 
Ve = 3000 km s"', and a=\ for all of their SNe, with no co- 
variance between any of these parameters. Because models 
with more ^^Ni need less gamma-ray trapping to produce the 
same bolometric luminosity at a given time, there is a large 
fit ting covariance betw een the ^^Ni mass and fo, mentioned 
in lScalzo et al\ (l2010h . The value of a is model-depend ent, 
and not a fundamental physical quantity, but as noted in 33.21 
above, a = 1 .2 (±0. 1) is also a common choice when no other 
prior is available from explosion models. Since a affects the 
nickel mass, a smaller assumed value of a results in a larger 
^^Ni mass, but a smaller ejected mass, as interpreted from a 
given bolometric light curve. In a self-consistent choice of 
parameters, and q will each depend in part on the mass of 
^^Ni and therefore on a. 

The value of Ve is difficult to measure directly, since ob- 
served velocities of absorption line minima may depend 
on temperature as much as density. Since appears 
squared in Equation [3] its contribution to the error bud- 
get on MwD is potentially quite large if treated as an in- 
dependent input. However, its value can be constrained 
within a range of ±300 km s"' by requiri ng energy conserva- 
tion. Following practice in the literature (iHowell et a/.ll2006l : 
iMaeda & Iwamotolllool . we calculate the kinetic energy Ek 
as the difference between the energy E^j released in nuclear 
burning and the gravitational binding energy Ec, and then set 

y, = (£^/6MwD)'/'. 

Calculating the energy budget of a SN la requires us to as- 
sume a composition. Our model considers four components 
to the ejecta: 

• ^^Ni, which contributes to the luminosity, E^, and Ec', 

• Stable Fe-peak elements ("Fe"), which contribute to Ef^ 
and Ec; 

• Intermediate-mass elements such as Mg, Si and S 
("Si"), which contribute to E^ and Ec', 
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• Unburned carbon and oxygen ("C/O"), which con- 
tribute only to Ec- 

The input parameters, which we vary using a Metropolis- 
Hastings Monte Carlo Markov chain, are the white dwarf 
mass MwD, the central density pc (needed in the calculation of 
Ec), the parameter a from Arnett's rule, the bolometric rise 
time frise.boh and the fractions /ni, /fb, and fsi of ^^Ni, stable 
Fe, and intermediate-mass elements within Mwd- We fix the 
fraction of unburned carbon an d oxygen /co = 1 - /pe ~ /m - 
/si. We use the prescription of iMaeda & Iwamotol (120091) to 
determine E^: 

[l.74/Fe+1.56/Ni+1.24/si] x lO^erg. (4) 

We use the binding energy formulae of lYoon & Langeil(l2005h 
for Eq = Ec{M^f^D,pc), where is the white dwarf central 
density. These ingredients determine Ve- We apply Gaus- 
sian priors a = 1 .3 ± 0. 1 (as for SN 2007if: IScalzo et a/.l2010l) 
and frise.boi = 20 ± 2 days, and on E{B- V)host according to the 
"sha llow TBC" val ues i n Table|4l We adopt = 0.025 cm~^ g 
after lJefferyl (fT99l and lStritzinger et ali (l2006h . 

One limitatio n with our approach is the use of Eq from 
lYoon & Langed (l2005h . derived for supermassive, differen- 
tially rotating white dwarfs. This formula remains an eas- 
ily accessible estimate in the literature for the binding en- 
ergy of a white dwar f over a wide range of masses, used by 
several other authors (iHqwell ef alll2006tljefferv ef a/.ll2006i: 
IMaeda & Iwamotol l2009l) ~ The models of lYoon & Langen 
(I2005h have bee n criticized on the grounds that they may not 
exist in nature d Pirol |200"8|). nor explode to produce SNe l a 
if they do exist (jSaio. & Nomotoll200llPfannes et a/.ll2010al) . 
However, it seems reasonable to assume that such models 
could represent a snapshot in time of a rapidly rotating con- 
figuration, such as that encountered in a white dwarf merger, 
which then detonates promptly rather than cont inuing to exist 
as a st able obj ect. The merger simulations of iPakmor et ali 
(12011!) and Pa kmor e t al. (2012), though they produce com- 
paratively little ^^Ni , show that prompt det onations in violent 
mergers can occur iPfannes et al.l (l2010bh simulated prompt 
detonations of rapidly rotating white dwarfs with masses up 
to 2. 1 Mq, and found that the amount of ^^Ni produced could 
be as high as 1.8 Mq, similar to SN 2007 if. 

In summary, our procedure directly extracts from a fit to the 
bolometric light curve only the quantities M^i and to (Equa- 
tion |2]i, and marginalizes, in effect, over a, fnse.boi, Ve, and 
q. The front end of our modeling technique varies the mass, 
composition, and structure of the SN progenitor as physical 
quantities which we wish to constrain. We then convert these 
inputs into physically motivated priors on the values of Vg and 
q, using Equations|2l[3] andlH and finally calculate the ejected 
mass Mwd- 

4.2. Including the Effects of a Shell 

The above considerations all apply to conventional 
exponential-equivalent models of expanding SN la ejecta. To 
explain the velocity plateaus of the SNe in our sample, how- 
ever, our model has a disturbed density structure where the 
high-velocity ejecta (included in the mass Mwd which under- 
goes nuclear burning) are compressed into a dense shell of 
mass Msh, traveling at velocity Vsh. In tamped detonation mod- 
els, such as the explosion models DET2ENV2, DET2ENV4 
and DET2ENV6 (Khokhlov et a/.lll993l:lHoflich & Khohklovl 



119961 hereafter "DET2ENVN"), such a shell is formed at the 
reverse shock of the interaction of the ejecta of an otherwise 
normal SN la with a compact (^ 10'° cm) envelope of mate- 
rial with mass Menv {external to, and not included in, Mwd)- 
The suffix in DET2ENVN refers to the envelope mass, so 
for example model DET2ENV2 has Me„v = 0.2 M©. "Pulsat- 
in g delayed detonation" models, such as the PDD535 model 
of iHoflich & Khohklovl (119961) . have similar shells created 
by non-homologous pulsations of the white dwarf progenitor 
prior to the final explosion, and hence do not require an exter- 
nal shell of material. However, these models tend to produce 
fainter events, with much shorter rise times and redder colors, 
than we observe for our sample, and so we do not consider 
them here. 

In a tamped detonation, the material which will form the 
shell imparts its momentum to the envelope, which in the 
DET2ENVN models acquires an average velocity of about 
1 .5\'sh. The interaction ends within about the first minute after 
explosion, and the shell then expands homologously with the 
other ejecta thereafter We observe Vsh directly as the plateau 
velocity, allowing us to constrain Msh and Menv, and, indi- 
rectly, the kinetic energy scale Vg of the ejecta. For a given 
value of Ve and a measured value of Vsh, and neglecting the 
binding energy of the en velope, conservatio n of momentum 
gives (for more detail see iScalzo et a/.ll20Tol) 



Me, 



Vsh \ Ve J \ Ve 



M 



WD- 



(5) 



where Q{a,x) = '^{a,x)/r(a) is the incomplete gamma func- 
tion. We calculate /env = Menv/MwD and f^h = Msh/M^^D 
by solving Equation |5] numerically. In a double-degenerate 
merger scenario, the total system mass Mtot = MwD+Menv is 
then equal to the initial mass of the two white dwarfs under- 
going the merger. 

For these calculations, we use only the velocity of the 
v(Si II) plateau measured from our spectroscopic time series. 
We do not model the duration of the plateau or the behavior 
of v(Si II) after the plateau phase ends. While the detailed 
evolution of v(Si II) undoubtedly contains useful information, 
reproducing it would require detailed calculations of synthetic 
spectra which are beyond the scope of this paper. However, as 
long as we have enough measurements of \'(Si II) to show that 
a given SN exhibits plateau behavior, we can reliably mea- 
sure I'sh without knowing the opacity of the material in the 
shell. We will compare our observations to previous numeri- 
cal models and observations of SNe la with interacting shells 
in ED 

4.3. Including the Effects of Central Density on Fe Yields 

In IScalzo et al\ ( 120101) . the stable iron fraction /pe was al- 
lowed to vary freely. In this situation, /pe and /si are nearly 
degenerate, since the contribution per unit mass of Fe to the 
nuclear energy released is only about 40% higher than that of 
Si. However, since Fe is produced by neutronization in the 
densest parts of the ejecta during the explosion, a high value 
of /pe can greatly reduce q because the formation of a large Fe 
core displaces ^^Ni to a higher average velocity and a lower 
optical depth. It therefore becomes important to constrain Fe 
production in any model in which we attempt to calculate q. 

The DET2ENVN explosion models, on which our models 
are loosely based, were intended to describe the detonation of 
a low-density white dwarf merger remnant of mass 1.2 Mq 
inside envelopes of varying mass. The central density pc in 
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these models is 4 x 10^ g cm ^, substantially lower than typi- 
cal central densities of ~ 10'' g cm"^ of deflagrations and de- 
layed de tonation models in the literature jNomoto et a/.ll984t 
iHoflich & Khohklov 1996t iKrueger et al.ll2010l) . These mod- 
els have no stable Fe cores immediately after explosion, with 
XsfiNi = 0.9 throughout the region where ""^Ni and Fe are pro- 
duced; 

IKrueger et al.l (120101) investigated the effects of central den- 
sity on ^^Ni yields in 3-D simulations of detonations of 
Chandrasekhar-mass white dwarfs, averaging over an en- 
semble of realizations for each value of pc- They found 
/Ni/(/Ni+/Fe) = 0.9 for p, = lO** g cm-^ decreasing by 0.047 
on average for each 10^ g cm""* increase of pc thereafter 

Our model already includes the effects of the central den- 
sity Pc on the bindin g energy Ec, of the whi te dwarf, through 
the fitting formula of lYoon & Langed (l2005h . Although Mwd 
may be super-Chandrasekhar in our models, the overall ex- 
tent of nuclear burning should depend on the density, not 
the mass, and so we may consider extrapolating those results 
here. Since the link between /pe and /ni is statistical rather 
than deterministic, we do not attempt to calculate a definite 
fraction for each model. Instead, we calculate the ratio of 
^^Ni to total iron -peak elements , 77 = fy j / (/ni + /pe), as well as 
the relation from IKrueger et all (120101) : 



T^Ki-io = min <^ 0.9,0.959-0.047 



W g cm- 



(6) 



enforcing a Gaussian prior 77 = ryKrio ±0.1. For most models 
calculated, this results in only a small fraction of stable Fe, as 
appropriate for low-density explosions. 

4.4. Calculation of q for SN la Models with a Shell 

The gamma-ray transport form factor q is the dimensionless 
^^Ni-weighted gamma-ray optical depth through the ejecta. 
Its value ranges between and 1, with large values corre- 
sponding to high central concentrations of ^^Ni; the gamma- 
ray optical depth is proportional to q. For perfectly mixed, 
exponentially distributed ejecta where the fraction of ^''Ni is 
cons tant throughout, q = 1 /3. 

In lScalzo et al.\ (120101) . the calculation of the total mass for 
SN 2007if assumes q = 1/3. We chose this value because 
the lack of a distinct second maximum in the /-band light 
curve suggested that a lar ge a mount of ^^Ni had to be mixed 
to higher velocities (see 33.1l above). This is not necessarily 
true for our other SNe. The use of ^ = 1/3 also assumes that 
the reverse-shock shell has a negligible effect on gamma-ray 
trapping, which may also not be true for very massive shells. 
Fortunately, q is easy to calculate numerically (Jeffery 1999); 



(7) 



/(, dzz^p{z)X56^Xz) 

where z and z' are dimensionless velocity coordinates in units 

of Vg, Zs = \/z'^ + z^ + 2zz'/i is the beam path length, p(z) is 
a dimensionless density profile normalized to unit mass, and 
^56Ni(z) is the (velocity-dependent) ^^Ni fraction just after ex- 
plosion. The geometry of the integration is shown in Figure|6] 
To calculate q, we assume a density profile motivated by 
DET2ENVN which includes an envelope with density profile 
p ^ z~^ on the outside, a thin Gaussian shell at Zsh = v'sh/v^, 
and an exponential in velocity beneath the shell. The total 




Figure 6. Integration geometry for the Ni distribution form factor 
q. A '^Co gamma-ray is emitted from point z in the region where 
^■561^1(2) > (inner shaded region), and travels along the ray Zs to- 
wards the observer. Point z! lies in a mass shell of thickness dz! 
(shaded annulus) in which the gamma-ray is trapped, scatters and 
deposits its energy. 



density is given by 

p{z) = Cbulk 

+ Cshexp 



Z-Zsh 



2 V fz,sh 



exp(-z) + Ce, 

21 



Zsh-Z \ _2 
Z 

CTz.sh 



(8) 



where we choose 0(z) = erfc(z/-\/2)/2 as a smooth function 
which approaches 1 for z ^ -1 and for z ^ +1, and the 
constants Cbuik, Cenv, Csh, are determined so that the mass 
fractions in eac h density profile c omponent agree with the 
input values. Cjui mby et al.\ (120071) suggest a nominal width 
of 500 km s"' for the reverse-shock shell at velocity Vsh ^ 
10000 km s~'. The se lf-similar shock interaction model of 
IChevalier et al.\ (119821) suggests that the reverse shock ve- 
locity width should be about 3% of the velocity at the con- 
tact discontinuity for an interaction with a envelope and 
ejecta with a power-law density profile r~" with n ^ 1 (used 
to approximate SNe la in the context of interaction w ith a 
CSM wind, e.g.. IWood-Vasev. Wang & Aldering||2004l) . We 
assume a shell half-w idth of a-^.sh = 0.015zsh in line with 
IChevalier et al.\ ( 119821) . although the expansion is homolo- 
gous and no longer self-similar after the interaction. Our re- 
sults are not sensitive to the exact value of CTz shi values in the 
range (0.01-0.05)zsh give us the same value of q to within 5% 
for reverse-shock shells with masses up to 0.5Mwd, and with 
much better agreement for less massive shells. 

For X56is,i(z), we use a parametrized composition structure 
inspired by Kasen (2006): 



X56Ni(z) ■■ 



( niFe.i 



- m(z) 



^Fexoie 

m(z)-OTiPE \ 

fliPE / 



X max(?7,0.9), 



(9) 
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Table 7 

Extracted model parameters 



SN Name 



M,o,/M0 (> 98%)" Mv/d/Mq (> 98%)'' M^-JMq" 



fen 



/Fe/sh « 98%)' X 



as 



SNF 20070528-003 
SNF 20070803-005 
SN 2007if 
SNF 20070912-000 
SNF 20080522-000 
SNF 20080723-012 



(> 1.49) 

-0.16 ' 

2.30;0;g (> 1.94) 



1.56+ 



(> 1.50) 
I (> 1.34) 



(> 1.32) 
1.38+« |5 (> 1.18) 
1.98+0-21 1 

(> 1.33) 
lA5t°-" (> 1.25) 



l.79^,fi (> 1.49) 1.69!j]-2| (> 1.41) 



0.91 ±0.12 
0.75 ±0.09 
1.37±0.13 
0.89±0.15 
0.80 ±0.11 
0.84 ±0.11 



0.08 ±0.05 
0.14 ±0.03 
0.16 ±0.04 
0.09 ± 0.05 
0.08 ±0.02 
0.06 ± 0.02 



0.30±0.12 
0.40 ±0.04 
0.44 ±0.05 
0.31±0.12 
0.28 ±0.04 
0.25 ±0.05 



0-23Lo.i4 
0.40!«:1^ 

+0.17 



+"■2'' « 0.66) 
(< 0.68) 



1.655™"(> 1.39) 
(> 1.77) 



1.50+5™1> 1.28) 
(> 1.46) 



-0.13 ' 

0.33^ \l (<0.59) 
« 0.45) 
0.30:"^;{^ «0.60) 
0.47+|;;|^ « 0.83) 



1.08 
0.69 
1.32 
0.10 
0.37 
1.09 



0.290 
0.654 
0.077 
0.775 
0.976 
0.324 



SN 1991T 
SN 2003fg 



0.77±0.15 
1.18±0.16 



0.10 ±0.03 
0.16 ±0.07 



0.34 ±0.05 
0.44 ±0.10 



0.83 0.600 
0.09 0.945 



Note. — Quantities with error bars are marginalized over all independent parameters. Uncertainties are 68% CL (Icr) and represent projections of 
the multi-dimensional PDF onto the derived quantities. Upper or lower limits on poorly constrained properties are 98% CL. 
" Total system mass of the two merged white dwarfs. 

l" Mass of the central white dwarf merger product which undergoes nuclear burning in the explosion. 
^*Ni mass synthesized in the explosion. 
Ratio of envelope mass to central merger product mass. 

Fraction of mass of burnt ejecta which is compressed into the reverse-shock shell. 
' Mass fraction of iron-peak elements ('*Ni -I- stable Fe) in the reverse-shock shell. 
^ Minimum xi achieved by a fit to the data. 

Probability of attaining the given value of or higher if the model is a good fit to the data, incorporating all priors. 
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Figure 7. Density and composition for an example model with /pe = 
0. 12, /n, = 0.6, /sh = 0. 12, /e„v = 0.02. Top: Mass of a shell of ejecta 
of thickness dz as a function of z. Black dotted line: overall density; 
solid lines, in rough order of increasing z: stable Fe-peak (blue), ^*Ni 
(red), IME -I- unbumed C/0 (green), C/0 envelope (black). Bottom: 
Composition of model as a function of mass coordinate m(z). 

where m{z) = /J z^piz) dz is the mass coordinate, and the stable 
Fe-peak core and the ^^Ni mixing zone are bounded by mass 
coordinates mpe.core and miPE with mixing widths flpe.coi-e and 
OjpE, respectively. This allows for some stable Fe-peak ele- 



ments to be mixed throughout while maintaining a core in the 
innermost regions for high central density models (see Fig- 
ure |7]i. For our modeling below we choose ape = 0.025 and 
flNi = 0.35, so that the outward mixing of ^^Ni c orresponds 
roughly to the "enhanced mixing" case of iKaseni (2006) , in 
accordance with the behavior we see in the light curves. As 
the mass of the shell increases, the results may depend more 
sensitively on the particular distribution of ^^Ni in the shell. 
We proceed with the analysis, but caution that more detailed 
models of the shock interaction, and/or actual hydrodynamic 
simulations of the explosion, may be needed to accurately 
understand gamma-ray transport for cases in which a large 
amount of ^^Ni is swept up into the shell. 

In general, the values of q are higher for our SNe (0.45 ± 
0.05) than the nominal q = 1/3 value for a completely mixed 
exponential SN la, but there is very little variation with shell 
mass fraction. The mixing of ^^Ni into a shell (and potentially 
above the photosphere) and the displacement of ^^Ni to higher 
velocities by a stable Fe core have comparable effects on q, 
but the former effect is minimized in the "enhanced mixing" 
model characteristic of our light curves. 

4.5. Modeling Results 

The results of our modeling for the SNe in our sample are 
summarized in Table]?] Figure]8]shows the confidence regions 
in the MwD-''''Ni mass plane for all six SNfactory SNe. While 
we limit our later statemen ts ab out relative rates (see 35.4l i 
and Hubble residuals (see ^5.61 ) to the untargeted SNfactory 
search, which samples the smooth Hubble flow (z > 0.03) and 
samples all host galaxy environments in an unbiased manner, 
we also include our results for two spectroscopically analo- 
gous SNe la from the literature as useful points of compari- 
son: SN 1991T and SN 2003fg (see Appendix]Ali. 

Our results can be summarized as follows for the SNe in the 
SNfactory sample: 



We recover the results of iScalzo et al.\ ( 120101) for 
SN 2007if, within the uncertainties. The total sys- 
tem mass has come down slightly, from 2.41 Mq to 
2.30 Mq (probability distribution median), since trap- 
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Figure 8. Confidence regions in the MwD- Ni mass plane for our six SNe, with priors a= 1 .3 ± 0. 1 and rj = 0.90 ± 0.05 as for DET2ENVN. 
SN 1991T and SN 2003fg are also shown. Contours bound regions of constant probability density. Colored regions are 68% (red), 90%, 95%, 
99%, and 99.7% CL (blue). 
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ping of Co gamma-rays by the envelope is now in- 
cluded in the mass estimate. The fractions of the total 
system mass in the shell and in the envelope, set by the 
plateau velocity, remain the same. 

• SNF 20070803-005 and SNF 20080522-000 are con- 
sistent with being Chandrasekhar-mass SNe la, albeit 
spectroscopically peculiar ones. 

• SNF 20070528-003 and SNF 20070912-000 have data 
quality sufficient only to obtain lower limits on the 
total mass, and upper limits on the shell and enve- 
lope mass. In particular the lack of late-time photom- 
etry points mean that the lower limits are driven by 
the ^''Ni mass from Arnett's rule. The total system 
mass Mtot, which includes the mass Menv of the enve- 
lope which we infer from the plateau velocity, is super- 
Chandrasekhar-mass at > 98% confidence, although 
the mass Mwd of the central merger remnant is con- 
sistent with Chandrasekhar-mass values. 

• SNF 20080723-012 appears to be a moderately super- 
Chandrasekhar-mass object, with ^^Ni mass, total mass, 
and other parameters intermediate between SN 2007if 
and the rest of the population. We place a 98% CL 
lower limit of 1.41 Mq on Mwd, so that the progenitor 
system is likely super-Chandrasekhar-mass even with- 
out including Menv 

We also perform the following cross-checks. First, our re- 
sults are not strongly sensitive to the assumption of a particu- 
lar mixing parameter, changing by less than 2% if we instead 
assume a completely stratified composition with ope = flNi = 0. 
Although changes in the mixing parameters do influence q 
substantially for models with small amounts of ^^Ni, they 
matter much less when the ^*'Ni mass is large. We also con- 
firm that when we discard the velocity information and fix 
Menv = 0, the resulting median reconstructed masses Mtot are 
within 2% of the median values of Mwd, and the probability 
distributions have comparable widths. This is expected, since 
only the most massive envelopes (as in SN 2007if) make an 
appreciable contribution to the gamma-ray optical depth as 
seen from the inner layers of ejecta. 

5. DISCUSSION 

We have shown that the subset of SNfactory events selected 
based on initial spectra similar to SN 2003fg, SN 2007if and 
SN 199 IT can be fit well by a simple semi-analytic model of a 
tamped detonation, intended to describe the results of mergers 
of double-degenerate systems with total mass at or exceeding 
the Chandrasekhar limit. In the following section we explore 
results in the literature which can help us determine the extent 
to wh ich t his interpretation is unique and relevant. 

In ^5.11 we compare our observations to numerical models 
of SNe la with interacting shells in their ejecta, and to ob- 
servations of SN 2005hj taken bv lOuimbv et al.\ (120071) . In 
^5.2[ we evaluate the possibility that some of the luminosity 
of ou r events may be due to an ongoing shock interaction. In 
^5.31 we discuss other recent models of asymmetric, single- 
degenerate SN la explosions which attempt to reproduce the 
velocity plateau phenomenon, and their implications both for 
whether our SN e are super-Chandrasekhar and whether they 
are mergers. In ^5.4[ we present the relative rates and examine 
the extent to which our results may constrain merger scenar- 
ios if a significant number of SNe la have double-degenerate 



progenitors. In ^5.51 we compare the aggregate probability 
distribution of total system mass from our events to popula- 
tion synthesis models of super-Chandrasekhar-mass double- 
degenerate mergers. Finally, in ^5.61 we examine the devia- 
tion from the Hubble diagram for these SNe and discuss im- 
plications for SN la cosmology. 

5.1. Comparison with Shell SN la Models in the Literature 

Our most recent point of comparison in the literature for 
the application of explosion models with in t eractin g shells 
to observations of SNe la is lOuimbv et al\ ( 120071) . They 
noted the presence of a velocity plateau in their observa- 
tions of SN 2005hj and compared them to delayed deto- 
nation, p ulsating delayed detonation and tamped detonation 
models (iKhokhlov ef a/.l [1991 iHoflich & Khohklovl [19961: 
iGerardv et a/.ll2004l) . 

For a shell mass fraction /sh = 0.07, similar to our less 
extreme SNe, the pr edicted B-V color is around 0.05-0.1. 
lOuimbv et al\ (120071) note that the systematic uncertainty in 
the absolute value of {B - y),nax may be as large as 0. 1 mag 
for the DET2ENVN models and other shell models. The color 
they report for SN 2005hj is (B-y)max = 0.04 ±0.06, consis- 
tent with the mean color in our sample, for a claimed shell 
mass fraction /s h = 0.14, comparable to SNF 20070803-005. 

According to lOuimbv et al.\ ( 120071) and references therein, 
one might expect to see cooler photospheres, and hence redder 
(B-y)niax, with increasing shell mass fraction. Our modeling 
predicts that SN 2007if, with the lowest plateau velocity and 
the most massive progenitor, has the most massive shell in 
both absolute and relative terms. This SN may be somewhat 
redder than the others near maximum light, but the lack of 
light curve coverage near maximum makes it difficult to say 
exactly how much. Table |7| also shows that SN 2007if has the 
largest expected fraction of iron-peak elements in its shell, 
corresponding to the material near the photosphere around 
maximum light. To the extent that SN 2007if is intrinsically 
redder than our other SNe, this may be due in part to line blan- 
keting by iron-peak elements, rather than a low-temperature 
photosphere, which would be inconsistent with the weakness 
of Si II A5800 in all the SNe in our sample. 

While we cannot predict from theory how long the plateau 
phase should last without more sophisticated modeling of our 
SNe, the durations of our plateaus are also broadly consistent 
with expectations. For models of Chandrasekhar-mass events 
with envelope masses about 0.1 MQ(Ouimby et al. 2003), the 
plateau phase is expected to last about 10 days. In most of 
our SNe it lasts at le ast 15 days (and at least 10 days for 
SNF 20070912-000). lOuimbv et al] dlOOTI) derive a plateau 
duration of 20 ± 10 days for SN 2005hj. 

SNF 20080522-000 shows the longest-Hved plateau in our 
data set, stretching as early as 10 days before B-band max- 
imum and lasting as long as 30 days. A more conservative 
estimate would be that the plateau phase is confirmed to last 
from day -3, when Si II A6355 becomes strong enough that 
the error bars on the velocity of the absorption minimum drop 
below 500 km s~', to day +15, where Fe II lines begin to de- 
velop near Si II A6355 and blending may become a concern. 
All of our velocity measurements in this 18-day time window 
are contained in a narrow range just 250 km s"' wide. 

5.2. Constraints on Ongoing Extended Emission 

iFrveref a/1 ( [2OT0l) ran three-dimensional smoothed-particle 
hydrodynamics (SPH) simulations of double-degenerate 
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mergers. They found that the central merger remnants are 
indeed surrounded by an envelope with an approximate ra- 
dial density profile p(r) oc with 3 < /? < 4. They then 
performed radiation hydrodynamics simulations of the inter- 
action of the SN ejecta with the envelope, calculating syn- 
thetic spectra and light curves of the resulting explosions. 
For sufficiently massive envelopes (more massive than about 
0.1 Mq), energy advected by the shock is released over the 
evolution of the SN, producing non-negligible luminosity near 
maximum light and extending emission into UV wavelengths. 
Based on these findings, J'ryer et al. (2010,) argued that these 
"enshrouded " systems would, in all likelih ood, look nothing 
like SNe la. iBlinnikov & Sorokina I (120101) performed analo- 
gous calculations using the STELLA radiation hydrodynamics 
code, finding that in general a radial density profile as steep 
as r~'^ looked similar to a normal SN la in optical wavelengths 
(UBVRI), whereas a profile varying as would be domi- 
nated by shock emission. 

These findings put significant constraints on the radial ex- 
tent and density profile of any envelope which might have en- 
shrouded the progenitors of the SNe in our sample. In partic- 
ular, the blue, mostly featureless spectrum seen in SN 2007if 
at phase -9 days is characteristic of what lFrver et al.\ (120101) 
expect, and the optical emission could be powered in part 
by advected heat energy from the shock interaction at early 
times, resulting in a broader light curve. However, the fact 
that v(Si II) in our events is seen to assume its plateau value as 
early as a week before maximum light, and the colors appear 
similar to those of normal SNe la near maximum light, argue 
strongly against any ongoing interaction with an extended en- 
velope or wind.^ 

As noted in IScalzo et al.\ (120101) , if some fraction of the 
maximum-light luminosity is due to shock heating, whether 
advected or resulting from a fresh interaction, instead of 
^^Ni decay, our mass estimates would tend to increase. 
This is because the influence of the shock interaction would 
diminish substantially by about 50 days after explosion 
(IBlinnikov & Sorokina 2010), and more ^^Co gamma rays 
would have to be trapped in order to reproduce the observed 
light curve. Any model in which shock heating produces a 
large amount of luminosity more than 50 days after explosion 
would probably not look like a SN la, and could not explain 
the appearance of the SNe la discussed in this paper 

One straightforward way to address the extent of any on- 
going consequences of shock interactions in future studies 
would be to obtain early-phase UV light curves of a candi- 
date super-Chandra SN la from a satellite such as Swift. The 
signature of any strong influence of shock heating would be 
evident therein. 

5.3. Implications of Possible Asymmetry 

The large inferred masses of our candidate super-Chandra 
SNe la, taken together with the observed plateaus in the Si II 
velocity, are both naturally explained by the tamped detona- 
tion model we put forth in this paper, and by the underlying 
double-degenerate merger scenario it represents. Some recent 
work, however, points towards the possibility of explaining 
these events in terms of asymmetric single-degenerate explo- 
si ons. 

iMaeda et al.\ (l2009l I2010allbl lIoTlh invoke large-angular- 
scale asymmetries to explain the diversity of velocity gradi- 
ents observed in normal SN la explosions and the low/high ve- 
locity gradient (LVG/HVG) dichotomy (.Benetti ef a/. 2005.) . 



They suggest that an asymmetric explosion may cause an 
overdensity in the ejecta on one side of the explosion, caus- 
ing LVG behavior when viewed from that side, with the other 
side relati vely less dense and e xhibiting HVG behavior Ad- 
ditionally, |Hachisu7raI](|20Tll) suggested that optically-thick 
winds blown from an accreting white dwarf could strip mass 
from the outer layers of its donor star, regulating the accre- 
tion rate and potentially allowing the white dwarf to accrete 
without exploding until reaching masses as large as 2.7 Mq. 
Such a white d warf would hav e to rotate differentially, as 
in the models of lYoon & Lange r (2005), with the attending 
uncertainty in the evolutionary history of such objects. If 
a differentially rotating white dwarf were to explode asym- 
metrically, this could present an explanation for our observa- 
tions within the single-degenerate scenario. A scenario like 
this one, if correct, could also explain SN 2007if and the 
HVG SN 2009dc as being similar objects viewed from dif- 
ferent angles. Tanaka et al. (2010) interpreted the low contin- 
uum polarization of SN 2009dc as evidence for a nearly spher- 
ical explosion, but low polarization could also be observed in 
an axisymmetric explosion viewed along the symmetry axis, 
as those authors note. 

While we cannot at this time conclusively rule out the pos- 
sibility that the Si II velocity plateaus we observe in our SNe 
result from asymmetry, no asy mmetry is needed a s yet to ex- 
plain them, as argued e.g. by iMaeda & Iwamotol ( 120091) for 
super-Chandrasekhar-mass explosions. The physical cause of 
the plateau — overdensity in the ejecta — is the same in sym- 
metric and asymmetric models. If the velocity plateaus are in- 
deed the result of asymmetric explosions, the shell mass frac- 
tions derived in Table |7] would still have meaning in terms of 
disturbances to the density structure along the line of sight, 
but the inferred envelopes would not be present, i.e., Menv 
would be zero and Mtot would equal Mwd for the SNe ana- 
lyzed in this paper. 

We can use the relative rate of supernovae spectroscopically 
similar to those in our sample (see ^5.4| below) to make some 
general statements about how well they can be explained by 
lopsided asymmetric explosions. If the SNe la in our sam- 
ple belonged to the same population as normal SNe la, and 
the spectroscopic peculiarity and brightness were due entirely 
to viewing angle effects (seeiJCasen 2004. for an asymmetric 
model for which this is true), a relative rate of ~ 2% would 
imply a range of viewing angles no more than ^15 degrees 
from the symmetry axis. For models in which the asymme- 
try does not tr anslate into a very peculiar spectrum, as seems 
likely for the iMae da efa Z. (2010b, 2011:) models with high 
^^Ni mass, we should look instead at the diversity of veloc- 
ity gradients among spectroscopic analogues. Since all of the 
(spectroscopically similar) SNe la in our sample have veloc- 
ity gradients at the slowly-evolving extreme of what the model 
of Maeda et al. (2010b) claims to produce, it seems plausible 
that, within the context of this class of models, the density 
enhancements in the ejecta of our SNe are roughly isotropic. 

Asymmetries resulting in overdensities in SN la ejecta 
along the line of sight could also affect gamma-ray trapping 
for that line of sight only . M ost of the trapping happens in 
low-velocity ejecta (see ^4. 5b . so the form factor q should 
not be strongly affected; we think it unlikely that q, and our 
ejected mass estimates, could vary by more than 10%. Our 
^^Ni mass estimate s from Arnett ' s rule should also be robust, 
since according to IMaeda et al.\ ( 1201 II) . the effects of asym- 
metry on peak brightness are least pronounced for the most 
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luminous SNe la. Based on these considerations and those of 
the preceding paragraph, the super-Chandrasekhar-mass sta- 
tus of SN 2007if and SNF 20080723-012 seems secure to us. 

A different class of asymmetries might arise in bipo- 
lar, rather than lopsided, explosions. An axially symmetric 
merger geometry, featuring a disk of disrupted white dwarf 
material in the equatorial plane, could in principle result in ve- 
locity plateaus and peculiar spectra when viewed edge-on, but 
a more normal SN la appearance wh en viewed pole-on. Thi s 
possibility was recently suggested bv lLivio &Pringld(l20TTh . 
who also suggested that SNe la from merger events might 
show increasing diversity in their velocity gradients with in- 
creasing mass of carbon. In this scenario, SN 2009dc could 
be an event like SN 2007if vi ewed pole-on, consi stent with its 
low continuum polarization ( iTanaka et a/.ll2010l) : SN 2007if 
should then correspond to the edge-on case, and should show 
substantial polarization. iLivio & Pringld (1201 ll) make no pre- 
dictions for relative rates in different SN la subgroups, which 
may depend on details of the hydrodynamic interaction of the 
SN ejecta with the disk. Nevertheless, an occluded region of 
ejecta narrow enough to explain their rare incidence would 
subtend only a small solid angle at the source, with only a 
small fraction of the ejecta decelerated. This would lead to 
broader, more complex line profiles than the ones we observe. 
It therefore once again seems reasonable that our SNe should 
be more or less spherically symmetric, with only moderate 
large-angular-scale asymmetries. Furthermore, any asymme- 
tries compatible with our observations should influence only 
the inferred envelope masses Menv, and should leave our esti- 
mates of MwD intact at the 10% level, as we noted above. 

Future observations, including nebular spectra of candi- 
date super-Chandra SNe la for direct comparison with the 
abovementioned studies, may allow us to make more specific 
statements about their asphericity. Spectropolarimetry, as for 
SN 2009dc (Tanaka et al. 2010), may also be helpful. 

5.4. Constraints on Merger Scenarios from Relative Rates 

Mass modeling shows that all six of our SNfactory- 
discovered SNe la have a high probability of having super- 
Chandrasekhar-mass progenitors. Because ours was a wide- 
field search and our spectroscopy screening was conducted 
impartially, that is without pre-selection based on colors, lu- 
minosity, lightcurve shape, host galaxy environment etc., we 
may estimate the relative fraction of such objects. 

A total of 400 spectroscopically-confirmed SNe la were 
discovered by the SNfactory. We did not spectroscopically 
screen candidates found to be fainter on the discovery image 
than on any pre-discovery detections that may have existed, 
as such candidates were assumed to be after maximum light. 
Furthermore, we did not perform spectroscopic screening in 
the cases where a host galaxy had a previously -known redshift 
beyond z = 0.08, though such cases were rare. 

We consider two subsamples — one that is purely flux- 
limited and another that is further limited to z < 0.08. For 
both subsamples we consider only SNe la discovered at or 
before maximum light, when the characteristic spectral fea- 
tures of SN 20 03fg, SN 2007if and SN 199 IT are clearest 
jLi et a/.ll201 lb . The flux-limit is represented by the survey 
efficiency as a function of magnitude rather than as a fixed 
magnitude. A subset redshift-limited to z < 0.08 is also con- 
sidered as a cross-check, since here details of the survey ef- 
ficiency become unimportant and the spectroscopic typing is 
very secure. 

Among the super-Chandra candidates, all survive the pre- 



maximum and flux-limit cuts, while SNF 20070528-003 and 
SNF 20070912-000 are eliminated by the redshift cut. The 
sample cuts applied to the overall sample leave 240 SNe la in 
the flux-limited sample and 141 in the redshift-limited sam- 
ple. Due to their enhanced brightness and generally longer 
rise times, the super-Chandra candidates are enhanced by fac- 
tors of 1.30 ±0.06 and 1.11 ±0.01 in each of the two sub- 
samples. These enhancements are calculated using the con- 
trol time in the SNfactory search, accounting for both appar- 
ent magnitudes and light-curve widths of these SNe relative 
to normal SNe la. The assigned uncertainties are illustrative 
of the effect of a generous ±0.25 mag shift in the detection 
efficiency curve. As expected, this affect is small and has es- 
sentially no impact for the volume-limited subset. Including 
these correction factors, the rates are 1.9;!;q5% 
for the flux-limited and redshift-limited subsamples, respec- 
tively. The difference in these relative rates is consistent with 
Poisson fluctuations, including accounting for the SNe la in 
common. 

These are the relative rates under the assumption that all 
of our candidates have super-Chandrasekhar-mass progen- 
itors, as their estimated masses under the tamped detona- 
tion model indicate. Alternatively, these can be taken as 
the rates for SNe la classified as like SN 199 IT by using 
SNID. Although our 1 .9-2.6% relative rate is somewhat lower 
than the "volume-lim ited" rate of 9.4^.1 j% (5-day cadence) 
from iLi et al.\ (1201 Ih . we note that it is consistent with the 
~ 1% "Ia-91T" rate of lSilv erman et al. (2012) from the same 
search. The differen ce between the two is related to how the 
events are classified. iLi et al.\ ( 1201 Ih use classifications from 
the lAU Circulars, rather than a homogeneous set of spec- 
tra over a given range of wavelengths and li ght-curve phases 
subclassified with a single me thod; in fact. | Silverman et al.\ 
(120 12h imply that some of the iLi et al.\ (1201 ll) subclassifica- 
tions may be photo metric (based on light-curve width) rather 
than spectroscopic. iLi et al.\ (1201 lb also treat SNe with type 
"la-99aa" as being "Ia-91T' rather than breaking them into 
separate su bclasses, which led to th eir much higher rate. Vs- 
ing SNID. [ Silverman et ai\ (12012b classify only one of the 
\Uetal.\ (|2011|) 1991T-likes (SN 2004bv) as "Ia-91T" and the 
rest as "Ia-99aa" or normal, leading to an updated "volume- 
limited" rate of 1 .4^g {] for the iLi et al.\ ( 1201 lb sample, con- 
sistent with our rate. As a point of comparison from an un- 
targeted search similar to SNfactory's, rather than a search 
targeting known galaxies, we also consider the SDSS spec- 
troscopic sample dOstman et al.WloHh . Seventy-eight of the 
141 SNe la in the sample of lOstman et al.\ ( 1201 lb have spec- 
tra with light-curve phases at or befo re maxirnum lig ht, proof 
agains t the 1991T-Uke "ag e bias" (iLi et al.\ IMTu) . Using 
SNID. IOstman et al.\ (1201 lb type only two of these 78 as "la- 
9 IT", leading to a relative rate of 3.4'!^fj%, again consis- 
tent with our own rate. We note that our rate estimates have 
considerably smaller uncertainty than these other 1991T-like 
SNe la rates. The combination of all three surveys, using our 
redshift-Hmited rates, gives a net 1991T-Uke rate of 2.3;|;o; 

The statistical confidence intervals from our models in- 
dicate the possibility that not all candidates are super- 
Chandrasekhar-mass. The probability density extending be- 
low the Chandrasekhar mass represents the equivalent of 0.2 
SNe la for both subsamples. If we correct for this excess prob- 
ability, the rates are lowered by factors of 5% and 3% for the 
flux-limited and redshift-limited subsamples. A more precise 
correction would depend on the true parent mass distribution 
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function, but it is evident from the above calculation that such 
details are well below the Poisson uncertainties. 

Finally, we consider the rates under the possibility that the 
shell model is not appropriate. In this case there would be no 
contribution from the shell mass, and the best mass estimate 
would be that of the central merger remnant, Mwd- In this cir- 
cumstance SNF 20070803-005 and SNF 20080522-000 could 
be considered more likely to have been Chandrasekhar mass 
events. The rates resulting from their removal would then be- 
come 1.5;!;q4% ^'^'I 1-0^o4%. While we view this situation as 
unlikely, for reasons discussed above in ^5.31 the rates remain 
dominated by Poisson uncertainties rather than modeling as- 
sumptions. 

These relative rates are very low. Therefore, if double- 
degenerate mergers are to contribute significantly to the total 
SN la rate, most such mergers would need to escape our spec- 
troscopic selection. Moreover, the SNfactory has not found 
overlumino us SNe la in addition to those presented here or 
SN 2005gj (lAldering et a/.ll2006h that would qualify for our 
redshift-limited sample. 

While our events are overluminous wi t h larg e ^^Ni masses, 
theoretical mergers jPakmoref a/.ll2011l 120 12l e.g.) thus far 
predict a range of smaller ^^Ni masses, resulting in fainter ex- 
plosions with lower photospheric temperatures. As discussed 
in ^52] above, a delay in the release of radiant energy from a 
shock interaction could also contribute to heating the photo- 
sphere at early times in events like SN 2007if, so that mergers 
with less massive envelopes would have more normal spectra. 

Merger events with less massive envelopes should also have 
less disturbed density structures, and hence may not display 
the low Si II velocity plateau behavior we observe in our sam- 
ple. The small envelope mass fractions implied for the ma- 
jority of merger events would then represent a constraint on 
either the dynamics of the merger in the case of a prompt 
explosion, or the post-merger evolution of the system in the 
case of a delayed explosion triggered by accretion onto, or by 
post-accretion spin-down of, the central merger remnant. Un- 
ambiguous detection of a density enhancement signature for 
less massive shells would be difficult. Construction of their 
mass distribution function would require a parent sample of 
SNe la all followed spectroscopically such that any merger 
candidates could be selected by their velocity plateau behav- 
ior rather than spectroscopic peculiarity. Moreover, such a 
study would require spectra at very early phases (-10 days 
or younger), since a less massive reverse-shock shell would 
become transparent at earlier phases, and might not subsume 
most of the intermediate-mass elements in the explosion. 

Determining the relative importance of any of these factors, 
and thus understanding the true relation of candidate super- 
Chandra SNe la to the general SN la population, will require 
further work and new spectroscopic data sets. However, we 
can say that a range of envelope masses are probably needed 
to explain the range of plateau velocities we observe in our 
data. 

5.5. Progenitor Mass Distribution in the Double-Degenerate 

Scenario 

Figure |9] presents the aggregate PDF of the masses of the 
SNe in our sample, each weighted inversely by their control 
time in the SNfactory search, accounting for both apparent 
magnitude and light-curve width. Its lower edge is roughly 
consistent with the Chandrasekhar mass, with a tail towards 
higher masses. 
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Figure 9. Total progenitor mass PDF for our sample (thic k solid line), 
compa red with vari ous population syn thesis models of iRuiter et al\ 
( I2009h as shown in lFrver ef al\ l l2010h : aA = 1.0 (thin solid line), 
a\ = 0.5 (dotted line), and 7 = 1.5 (dashed line). Population synthe- 
sis models have been cut off below 1.4 Mq and convolved with the 
mass resolution of our reconstruction technique, then normalized to 
match our total event rate. 

This distribution is qualitatively similar to expectations 
for super-Chandrasekhar-mass merging white dwarf systern 
masses from population synthesis models (iRuiter et a/.ll2009l: 
Hryererfl/. 2010). The limited statistics and theoretical un- 
certainties make it difficult to draw a firm, quantitative cor- 
respondence between observations and theory, but we can 
nevertheless look for gross inconsistencies. Figure |9] there- 
fore also includes three theoretical distributions of the pro- 
genitor system mass (Mtot) shown in iFrver et al] (120 1 Ol) . cut 
off below 1 .4 Mq and convolved with a two-sided Gaussian 
{<T- = 0.08 Mtot, = O.lOMtot) to model our reconstruction 
precision. 

Our observations correspond well with the IRuiter et'al] 
(120091) models on the leading edge, particularly the 7=1.5 
model. The differences are more pronounced on the high- 
mass tail; the agreement may improve with a larger sample of 
SNe. The other two models agree less well, predicting more 
high-mass mergers than we see in our sample. In the 7=1.5 
model, the relative rate of events dN /dMiot is about a factor of 
3 higher at 1 .4 Mq than at 2. 1 Mq, consistent with our obser- 
vations. The other two models predict more high-mass merg- 
ers than we observe; we see no evidence for a sep arate forma- 
tion ch annel or peak near 2 Mq as suggested by iFryer ef all 
(l2OT0l) . If such massive mergers do exist, they must therefore 
produce less ^^Ni, and be less luminous at maximum light, 
than the SNe in our sample. 

5.6. Hubble Residuals and Implications for Cosmology 

Table [3] also contains the Hubble residuals of each SN in 
our sample from a standard ACDM cosmology: 

A/i = ;7Zb - 25 - 5 log 1 (£//. /Mpc) - a(s-\) + I3c-Mb, (10) 

where nig, s and c are deri ved from the S ALT2 fits, di is calcu- 
lated using lWrightl (120061), and we use val ues of a, (3, and Mb 
from the w = -1 fit of lSulIivan et alUlOl Ih . (WeuseMB=M^ 
i.e., the absolute magnitude appropriate for SNe in host galax- 
ies with stellar mass less than 10^° Mq.) 

Figure [To] plots the relation A/i vs. Mwd (instead of M^t, 
so that it does not depend on interpreting the velocity plateau 
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Figure 10. Hubble residuals vs. central merger remnant mass for our 
sample (black points). SN 1991T (red) and SN 2003fg (green) are 
also plotted for comparison. Error bars are 68% (la) uncertainties. 
The dotted black line is the best-fit linear trend. 

in terms of an envelope) for the six SNfactory SNe la. For 
comparison we also plot SN 199 IT and SN 2003fg, mod- 
eled in a similar fashion (see Appendix |A]i. The Hubble 
residuals show a correlation with increasing mass (Pearson 
r = -0.76, p = 0.03). SN 2007if is 1.3 mag overluminous 
for its stretch and color, and by 1 mag even if no correction 
is made for its redder-than-average color. SNF 20080723- 
012 is about 0.6 mag too luminous for its stretch and color, 
putting it off the Hubble diagram but to a lesser degree 
than SN 2007if. At the low-mass end, SNF 20070803-005 
and SNF 20080522-000, while slightly overluminous, are 
not alarmingly discrepant, being within two standard devi- 
ations of the ACDM Hubble diagram assuming an intrin- 
sic dispersion of 0.15 mag. The higher-redshift (z ^ 0.12) 
SNF 20070528-003 and SNF 20070912-000 are too luminous 
for their stretch and color by about 0.5 mag; the residual de- 
creases to about 0.4 mag if no correction is made for the color. 

If the correlation shown in Figure [TO] is real, and the large 
Hubble residuals of SN 2003fg and SN 2007if are related 
to their large progenitor masses, the large Hubble residuals 
of SNF 20070528-003 and SNF 20070912-000 provide fur- 
ther circumstantial evidence that they too may have super- 
Chandrasekhar-mass progenitors. One might expect devi- 
ations of super-Chandra candidate SNe la from the width- 
luminosity relation for normal SNe la, either because the 
large ejected mass could increase the diffusion time of radia- 
tion through the ejecta relative to Chandrasekhar-mass explo- 
sions, or because the enhanced mixing of ^^Ni into the outer 
layers of ejecta in th ese events affects the ion ization state at 
the photosphere (see iKasen & Wooslevll2007l and references 
therein). 

Previous authors have warned about overlum inous SNe la 
that could skew the cosmological parameters. iRei ndl et al.\ 
(l2005i) found that the SN 1991T subclass give a mean Hub- 
ble residual of -0.4 mag. The sample presented here has a 
mean Hubble residual is -0.56 mag, dropping to -0.41 mag if 
SN 2007if were excluded on the grounds that it would be easy 
to eliminate. However, without a spectroscopic veto at high 
redshift, or by relying only on photometry at any redshift, it 
will be difficult to tag and remove SNe such as these from 
the Hubble diagram. Although SN 2007if is excessively over- 
luminous and would probably be excluded from any Hubble 
diagram on that basis, somewhat less luminous analogues like 
SNF 20080723-012 probably could not be safely removed 



based on their overluminosity alone. It would also be diffi- 
cult, without spectroscopy, to distinguish intrinsically overlu- 
minous SNe la from normal SNe la magnified by gravitational 
lensing; since len sing conserves photons in general relativity 
(lWeinberg|II976l) . all events must be included on the Hubble 
diagram to avoid bias in the reconstructed cosmological pa- 
rameters. 

Suppose candidate super-Chandra SNe la occur at the rel- 
ative rate we calculate for our redshift-limited subsample 
in ^5.41 (2.6"'"() ^) and have a mean Hubble diagram residual 
of -0.41 mag, the mean of the Table [3] entries excluding 
SN 2007if. The mean absolute magnitude of SNe la is then 
0.01-0.02 mag higher (68% CL) than it would be if these SNe 
followed the normal relations. This numb er is already com- 
parable to the 2% error budget suggested in lKim et al.\ (l2004h 
for Stage IV SN la cosmology experiments, and would be re- 
alized if the rate evolves by either doubling or going to zero 
at high redshift. 

6. CONCLUSIONS 

We have observed a sample of SNe la which we selected 
on the basis of their pre-maximum similarity to SN 2007if, 
searching for new candidate super-Chandra SNe la. Based on 
these spectroscopic analogues, we find that SN 2007if proba- 
bly lies on the extreme luminous end, with a hot, highly ion- 
ized photosphere and a disturbed density structure resulting 
from a stalled shock within the ejecta. We interpret these char- 
acteristics as evidence for an interaction with a circumstellar 
envelope early in the SN's evolution, as one might expect in a 
double-degenerate merger scenario. 

To reconstruct the masses of the progenitors, we apply an 
update d version of the modeling technique of IScalzo et al.\ 
dioiOl) . By applying this technique to SN 199IT for com- 
parison with lStritzinger et'al] (120061) . we show that including 
NIR flux near phase +40 days as part of the bolometric flux 
has a significant impact on the final reconstructed mass; while 
there is some uncertainty associated with correcting for NIR 
flux not actually observed, ignoring such a correction causes 
the mass to be systematically underestimated. Our modeling 
technique includes a set of Bayesian priors motivated by the- 
oretical models, capturing the covariance between different 
parameters of the system, and we sample the full probabil- 
ity distribution using a Monte Carlo Markov chain. From our 
sample, SNF 20080723-012 now joins SN 2007if in having a 
super-Chandrasekhar-mass progenitor at high statistical con- 
fidence. The mean error-weighted system mass for the four 
remaining SNe in our sample (1.63 Mq ±0.10 Mq) is also 
above the Chandrasekhar mass at a statistical confidence level 
exceeding 99%. 

Although the spectrophotometric observations discussed 
here are quite detailed, with excellent temporal and wave- 
length coverage, the modeling presented is relatively simple 
and is meant to explore a large parameter space of rudimen- 
tary explosion models quickly and efficiently. Detailed com- 
parison of the spectra to synthetic spectra of a range of con- 
temporary explosion models will be necessary to determine 
the extent to which our modeling faithfully represents the den- 
sity structure found in these SNe. Our simple model also as- 
sumes spherical symmetry of the SN ejecta, and this approxi- 
mation may not hold true in a significantly aspherical merger 
scenario. The shell mass estimates must necessarily be sub- 
ject to systematic uncertainties of this sort, and hence might 
be more properly construed as upper limits for which equality 
holds in the case of spherical symmetry. The total mass esti- 
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mates, which concern the bulk ejecta and not only the outer 
layers, should be more robust. 

Further study of a larger sample of such SNe la may provide 
additional insight into both the explosions themselves and the 
formation of their progenitor systems, particularly the extent 
to which velocity plateaus are linked to super-Chandrasekhar- 
mass progenitors. Future data sets including not only optical- 
wavelength observations, but infrared observations and neb- 
ular spectra, can help to verify the luminosity, ^^Ni content 
and distribution, and total mass of these events. Ultraviolet 
and X-ray observations of SNe la spectroscopically similar to 
SN 199 IT at early phase may also help to constrain the con- 
tribution of prompt shock emission to their luminosity, and 
the density and spatial extent of any envelope of material that 
might surround their progenitors just prior to explosion. Po- 
larimetry and spectropolarimetry, as well as velocity offsets 
measured from nebular spectra, can help constrain the degree 
of large-angular-scale asymmetry. Finally, similar analysis 
of a much larger SN la data set may help to discover super- 
Chandrasekhar-mass SNe la whose spectroscopic appearance 
is unlike those in our sample; this may help to clarify the con- 
nection of these events to the general SN la population, with 
dividends paid to the understanding of SN la progenitors and 
to the precision of SN la cosmology. 
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APPENDIX 

MASS MODELING OF SPECTRAL ARCHETYPES 

In addition to SN 2007if, SN 2003fg and SN 1991T are spectroscopic archetypes for the candidate super-Chandra SNe la 
presented in this paper Thus, in order to have self-consistent mass estimates for the full archetype set and provide the reader with 
familiar points for comparison, we have applied our modeling procedure to SN 199T and SN 2003fg. 

SN1991T 

SN 1991T is the most well-known spectroscopic archetype of our sample, although not the most spectroscopically extreme 
example. SN 1991T has good wavelength coverage in published spectra taken between maximum light and phase +40 days, 
allowing us to use the same analysis tec hniques as for the rest o f our sample. In addition, modeling of SN 1991T allows a useful 
direct comparison to the earlier work of IStritzinger et a{\ (120061) . 

SN 199 it's remarkably flat velocity evolution was pointed out bv lPhillips et al\ (119921) . and it has the second-lowest velocity 
gradient v in the original sample of Benetti et al. (2005). The DET2ENV2 model fits SN 1991T (lH5flich & Khohklov 1996), 
and in deed SN 1991T wa s once suspected to have had a super- Chandrasekhar-mass progenitor, mostly becau se of its high lumi- 
nosity (iFisher et a/.lll999l) . Later measurements of the distance (iRichtler et aZ.ll20()lt iGibson & Stetsonll2001h to the host galaxy, 
NGC 4527, have decreased the SN's luminosity considerably, and hence its inferred mass. 

We can make a new estimate of the progenitor mass of SN 199 IT using our reconstruction mefliod, which considers the 
gamma-r ay trapping in more detail than previous studies. We use the BVRI light curve of iLiral (119981) and the spectra l time 
series of iMazzali etal.\ ( 119951) obtained through the SUSPECT online supernova spectrum archive~ |Kanbur TTai] (12003) give 
the Cepheid-based distance modulus to NGC 4527, after correction for metallicity, as 30.71 ±0.13 (LMC period-luminosity 
relation) and 30.78 ±0.1 3 (Milky Way period-luminosity relation). We adopt the mean of these two measurements (30.74 ±0.1 3) 
for our work. From a joint fit to the spectral time series we extract £'W(Na I D)= 0.84 ±() . 1 8 A, givin g E(B- V)host ~ ■ ^ 3 ± 0.03 
from the shallow TBC relation, in good agreement with the estimates of iPhiUips et al.\ (119921) and iPhiUips et al.\ (119991) . Our 
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model represents a good fit to the available data, with = 0.83 (Pfit = 0.60). We recover a ^^Ni mass of 0.77 ± 0.15 Mq and a 
white dwarf mass 1.50;!;q [3 Mq. If we suppose that the plateau is caused by interaction with a uniform spherical carbon/oxygen 
envelope, this envelope has a mass of about 0.1 Mwd, and the total system mass estimate goes up to 1.65;!;o ^o- The results 
are similar to those for the nominal Chandrasekhar-mass SNe in our sample, SNF 20070803-005 and SNF 20080522-000. 

Stritzinger et al. (2006) measured a mass of 1.21 ±O.36M0 for SN 199 IT. Our median derived mass for SN 199 IT exceeds that 
of iStritzinger et al.i t2006i) bv over 0.4 Mq; our 98% CL lower limit onMwo is 1.28 Mq. However, we find that we can reproduce 
IStritzinger et al.\ (l2006l) 's numbers for fo, ^^Ni mass, and ejected mass if we neglect NIR corrections and ignore covariances in 
our fitting procedure, particularly that between q and v^. In particular, by including NIR corrections we find fo = 48 days, much 
larger than the[Stritzinger et al. (2006) value of 34 days and suggesting more massive ejecta. Our larger median value for q (0.42) 
and slightly lower Vg (2850 km s~') pull in the other direction, resulting in a mass closer to the the Chandrasekhar mass. 

In summary, SN 199 IT shows the observational hallmarks of a tamped detonation and may well have been a double-degenerate 
merger. Because of uncertainties in the host distance and reddening correction, we are unable to establish with confidence whether 
SN 1991T was itself super-Chandrasekhar-mass. 

SN2003fg 

SN 2003fg dHowell et a/.ll2006l) has a comparatively limited data set for our purposes, with only one spectrum at maximum 
light and no photometry at sufficiently late times for us to constrain the gamma-ray transparency. Nevertheless, there are noted 
similarities to SN 2007if: The single spectrum available is a good match to SN 2007if at a comparable phase, though with stronger 
Si II and S II. The observed /'-band light curve of SN 2003fg shows no significant inflection, like SN 2007if and unlike the less 
massive SNfactory SNe la. The low (^ 8000 km s~') Si II velocity near maximum light is also like SN 2007if at similar phases, 
and makes SN 2003fg a good candidate for a tamped detonation, although SN 2009dc showed a similarly low velocity without a 
discernable plateau ( Yamanaka et al. 2009). 

We can obtain a lower limit on the mass of SN 2003fg assuming it can be fit by a tamped detonation. We use the value of the 
bolometric absolute magnitude at maximum light (Mboi = -19.87) and the velocity of the Si II A6355 absorption minimum (~ 
8000 ±500 km s~') derived bv lHowell et al.\ (120061) . Assuming no host galaxy reddening, we find a 5^Ni mass of 1.18±O.16M0 
and a 98% CL lower limit on the total s ystem mass of 1 .IIMq (1 .46 Mq if the envelope mass is neglected). These are somewhat 
looser constraints than those of How ell et al\ ( 120061) . but are free of the assumption that the photospheric velocity should be 
related directly to the kinetic energy released in the explosion, which, as we have seen, may not be a good approximation for 
some density structures. 
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